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ABSTRACT 

The recent availability of timely synoptic earth 
imagery from the Earth Resources Technology Satellites 
(ERTS) provides a wealth of information for the monitoring 
and management of vital natural resources. The quantity 
and rate of acquisition of this data requires that the 
resource planning function be aided by computer processing, 
however, computer software systems have not existed to 
provide the planning function with an adequate computer 
information structure to maintain image interpretations. 
Available systems also lacked the simplistic interactive 
interface required by the resource planner who is not 
versed in the complexities of computer processing. 

Formal language definitions and syntax interpretaion 
algorithms have been adapted to provide a flexible, compu- 
ter information system for the maintenance of resource 
interpretation of imagery. These techniques are incorpora- 
ted, together with image analysis functions, into an 
interactive Resource information Management and Analysis 
System, IRIMAS, which is implemented on a Texas Instruments 
980A minicomputer system augmented with a dynamic color 
display for image presentation. A demonstration of system 
usage and recommendations for further system development 


are also included. 
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CHAPTER I 
INTRODUCTION 

Humanity is consuming ever greater quantities of land 
and natural resources; a situation which is causing great 
concern among many groups within our society. In this 
"land of plenty", shortages in desired commodities have 
already become apparent and the efficiency of the conver- 
sion of resources to the "good life", and even the "good 
life" itself, are being questioned and criticized. Modi- 
fications in the nominal mode of living in the United 
States, as in many parts of the world, are being dictated 
by resource shortages, public opinion, and legislative 
reaction. A major emphasis is being placed on the develop- 
ment and utilization of resource management systems to pro- 
vide for more efficient conversion of valuable resources, 
improved agricultural production, and minimization of the 
effects of blight on the environmental and social systems. 

A key element in the development of resource manage- 
ment systems is the acquisition and utilization of large 
quantities of information which is varied in nature and 
areal extent. Consequently, Earth Observations has emerged 
as an important field in the support of resource management. 
The techniques of acquiring the regional information have 

This dissertation follows the style of the Proceedings 
of the IEEE. 
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developed to the point that resource managers are now being 
overwhelmed by the great quantities of raw data generated 
by current technology. As a result, it is important that 
the technical knowledge developed by practitioners of the 
art and science of remote sensing be distilled into an 
analysis/display system amenable to utilization by resource 
management agencies at local and regional levels. It is 
only in this way that the vital information contained in a 
synoptic and timely overview of our planet can be applied 
immediately and specifically to the problems facing our 
society as a whole. 

It is the objective of this research to develop an 
information management and analysis system suited to the 
general needs of resource managers. The information system 
described implements the analysis techniques of the Earth 
Observations community in a minicomputer software/hardware 
system developed to provide a convenient, interactive 
interface between earth observations data and its users, 
and to implement a flexible multipurpose data handling 
scheme suited to handle the spatial and temporal informa- 
tion required for resource -related decision making. 

An Information System 

Several farsighted authors [1], and [2] have described 
the salient features of a computer land-use information 
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system. Such an information system could as easily apply 
to a broad range of resource management information require- 
ments. The foremost requirement of such an information 
system must be the effective and efficient presentation of 
data from multiple sources in the form of products which 
are relevant to the problems and processes of resource 
management . 

According to Brooner and Nichols [2], 

If any method is to be successful, it must 
meet the criteria of timeliness, accuracy, 
and economic feasibility. Data must be 
recorded, massaged, and retrieved in a 
manner that local and regional agencies 
that most desperately require such infor- 
mation can effectively utilize it. 

Condensing the thoughts of these authors, resource manage- 
ment information systems should consist of the following 
components: (1) system management, (2) data files, 

(3) interpretation/analys is , and (4) output product. 

The system management component provides a house- 
keeping/input function. This function provides for the 
input of data and its transformation into formats suitable 
for insertion into data files. It translates data from 
remote sensing imagery, surveys, and administrative records 
into machine records. Additionally, the system management 
component includes the software monitor of the total system, 
and constitutes the input interface to the resource 
management activity. 
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The data file component separates data into disaggre- 
gate data, areal-unit summaries, and time-period summaries. 
The disaggregate data, which includes satellite and aerial 
imagery, survey information, etc., is temporarily stored 
information related to a common geographic base. Areal- 
unit summaries, and time-period summaries include geograph- 
ically referenced resource information and other generally 
related parameter information which is further analyzed and 
displayed in a spatial format suitable for individual 
interpretation. These information files also provide the 
necessary data for the compilation of area statistical 
summaries . 

The interpretation/analysis component functions are 
not fully defined by Dueker and Horton [1], and Brooner and 
Nichols [2], but both sets of authors have recognized the 
ultimate necessity of these general functions in a resource 
management system. Swanlund [3], and Viglione [4], have 
described some of the current techniques utilized in inter- 
preting and handling image and spatial data. Numerous 
other authors [S], [6], [7], have described the theory and 
techniques of image enhancement and thematic classification 
from image survey systems. 

Additional authors have described statistical model 
analysis procedures for the interpretation and prediction 
of urban change and urban transportation requirements, but 
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many of these studies are limited in scope and interpre- 
tability because of the immense information management 
problem. Harris [11] and Candeub [12], have indicated that 
the development of computer information management systems 
and more adequate display of spatial information hold the 
key to further development in urban modeling and its use 
in the planning function. Many systems engineering study 
groups have recently become active in the application of 
systems analysis techniques to social and resource problems 
[13] > [1^3 • Yet again, the key to meaningful results is 

the management of the spatial information required by the 
models. 

One of the most significant aspects of an information 
management system is its output product, for it is upon 

this component that the whole system must be judged. With 
these words , 

In terms of land resource management 
the time-space aspect of the output 
product cannot be overemphasized. 

Ecology also requires graphic and 
cartographic methods especially 
suited to analysis of spatial data. 

Brooner and Nichols [IS] emphasize the output requirements 

of thematic maps and statistical displays, and Thomlinson 

[16] says that these type products are indispensable for 

ecological analysis. 

The works of many authors have served to provide a 
strong conceptual definition of a computer resource 
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information management system, and have provided the inputs 
from practitioners of the science of urban land management 
required for the engineering of such a system. Many other 
authors have provided the techniques required for the 
automatic and manual analysis of spatial and temporal 
inputs to a resource information management system afforded 
by the science of earth observations and remote sensing. 
Recent hardware developments in minicomputer technology and 
low-cost display systems have essentially provided the raw 
materials required for the development of a computer 
information management and analysis system amenable for 
utilization by a resource management function. 

Ob j ective 

It is the objective of this research to design and 
develop an Interactive Resource Information Management and 
Analysis System (IRIMAS), configured for simplistic user 
interaction and cost-effective implementation. This system 
will facilitate integration of remote sensor surveys and 
will provide access to this composite data base through 
spatial display media. Techniques for information enhance- 
ment and class identification will be implemented and 
included in the system to provide a versatile management 


tool . 
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Scope of Report 

The remote sensing field is varied and covers broad 
interdisciplinary areas. The specific equipment, tech- 
niques and physical considerations with regard to energy 
mechanisms which are being measured; the analytical and 
automated methods of analysis of data; and the application 
of results to diverse activities all tend to discourage a 
cohesive understanding of the techniques and philosophies 
involved. Chapter II is written to give a general overview 
of the field of remote sensing and remote sensing activi- 
ties covering energy mechanisms and spectral distributions, 
sensor systems in operational use, and philosophy of remote 
sensing analyses. It also suggests the needs of resource 
managers in the application of remote sensing techniques 
to problem solving. 

A resource information management and analysis system 
is described in a general manner in Chapter III. Format 
difficulties with current analysis systems and desirable 
characteristics of a comprehensive information system are 
described, and a modified format for resource information 
entry and retrieval are defined. Chapters IV, V, and VI 
discuss in technical detail the Interactive Resource 
Information Management System (IRIMAS) . A demonstration 
of the Interactive Resource Information Management System 
utilizing data available from Brazos County, Texas, is 
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discussed in Chapter VII, while conclusions and recommen- 
dations for system improvements are covered in Chapter VIII. 
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CHAPTER II 
BACKGROUND 

Remote Sensing 

Remote Sensing is an aggregate of activities concerned 
with inferring the nature of objects or areas by their 
measurement at a distance. These activities embody the 
physical sciences concerned with energy mechanisms, since 
it is generally changes in energy which are measured; the 
engineering sciences concerned with the detection and 
recording of the energy changes; and the analytical sci- 
ences, including the activities of mathematics, statistics, 
information and computing sciences; photogrammetry and 
image interpretation; and the other sciences which trans- 
form the recorded measurements. Also included are the 
application disciplines such as agriculture, geology and 
sociology which combine discipline experience with the 
information provided by the other groups to make the 
required inferences. The remote sensing field has evolved 
from a forced intermingling of these diverse activity 
groups; each with its own terminology, techniques, phil- 
osophies, and performance evaluation criteria. The worth 
of the composite remote sensing activity is judged by the 
social and economic impact of its products. 

Remote sensing activities have been modeled by Rouse 
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and Schell [17] in terras of a large scale information 
system directed toward the provision of specific applica- 
tion products. This system can be viewed in terms of the 
four functional areas represented: Energy Interactions, 

Energy Sensing Systems, Analysis Activities, and Management 
Evaluation and Action (See Figure II-l). 

The remote measurement of environmental characteris- 
tics is directly related to the measurement of energy 
surrogates of those characteristics. Energy released by 
a variety of mechanisms is interacted with by the 
Environmental Situation in a characteristic fashion. To 
the extent that these energy interactions are unique for 
a particular environmental situation, a knowledge of the 
character of the energy before the interaction and a 
measurement of its characteristics after the interaction 
can lead to a knowledge and characterization of the 
situation which produced the changes. 

Well developed remote sensor systems are available for 
the routine measurement of energy. Some of these sensor 
systems have been applied to a general data gathering 
mission and have provided a variety of measurements of 
the environment in a timely and systematic manner. 

Similarly, data handling and analysis techniques have 
been developed along well specified patterns and form 
identifiable groupings. It is primarily in the provision 



active 



Figure 1 1 - 1 . Earth Observations Information System 
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of data products to the Management Evaluation and Action 
activity and the specific products useful to an environmen- 
tal manager that the system is ambiguous. As such, the 
output data products produced by the analysis activities 
tend to be self-serving rather than fully beneficial to the 
next function block in the earth observations system. 

Energy Mechanisms 

The primary source of information available in remote 
sensor measurements is the spectral distribution of re- 
flected or emitted energy. Some comprehension of the 
sources and change characteristics of energy is necessary 
to obtain a fuller understanding of the measurements of 
remote sensing and their applicability in specific problem 
solving exercises. The full spectral range of energy from 
light well into the microwave region is measured by remote 
sensors, however, the majority of resource characterization 
by remote sensing techniques has been accomplished with the 
help of visible and near-visible radiation. It is this 
spectral region which is receiving primary emphasis, 
although microwave measurements and interpretation offer 
great future promise. 

Fundamental to an appreciation of energy changing 
mechanisms in the visible and near-visible regions is a 
knowledge of solar radiation and the atmospheric effects of 
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absorption and scattering. The characteristics of vegeta- 
tion, soils, and other natural resources which change the 
energy spectral distribution are important to the success- 
ful application of automatic classification algorithms, 
and also to manual interpretation of sensor measurements. 
The total energy change is characterized by a spectral 
signature which is the measured energy distribution influ- 
enced by the sensor target. Within the uniqueness of the 
characteristics of this signature often lies the success 
or failure of the measurement interpretation task. 

Electromagnetic Radiation and Spectral Color 

Two common models are utilized to describe electro- 
magnetic radiation, the wave model and the particle model. 
According to electromagnetic wave theory, energy is 
propagated through space in a wave pattern analogous to 
sound waves in the air or shock waves in water. An 
electromagnetic wave consists of two fluctuating fields 
which are generated as a result of an electrical charge 
being accelerated and the wave frequency depends on the 
number of accelerations per second. The second approach 
to the study of radiative energy transfer is explained by 
particle or photon theory. This theory states, in effect, 
that energy is transferred in discrete units called quanta 
or photons and that the energy of a photon varies directly 
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with frequency. Both theories are consistent and serve to 
provide a complementary understanding of electromagnetic 
energy radiation over the widest spectral range. 

The term "spectrum’' has been applied to the continuous 
series of colors obtained when white light is separated 
into its component colors of red, orange, yellow, green, 
blue, indigo and violet. This series taken together is 
called the visible spectrum. On the basis of the wave 
theory of light, color is associated with the radiation 
wavelength. In the visible spectrum, red has the longest 
wavelength and violet has the shortest wavelength. The 
visible spectrum is only a small portion of a much broader 
band of radiation exhibiting properties similar to those of 
visible light. This total range of radiation is called 
the electromagnetic spectrum. This enlarged spectrum 
definition includes radiowaves, microwaves, infrared or 
heat rays, visible light, ultraviolet rays. X-rays, and 
Gamma rays, although in remote sensing the utilization of 
microwave, infrared, visible and ultraviolet radiation is 
more significant. 

Variations in color represent the interpretation that 
man's eyes and brain make of the differences in wavelength 
or frequency of light. Color is always associated with 
the way things are seen. Two simple concepts help define 
the color of an object. If various colors of light are 
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allowed to fall on a white screen, the result is that each 
color is faithfully reflected, indicating that a white 
object reflects all colors or wavelengths. If white light 
is allowed to fall on a colored surface, say red, the 
surface absorbs all the colors except red, which is reflec- 
ted. Therefore, the color of a material can be defined 
as that of the light which it reflects. 

Absorptance is that quality or characteristic of an 
object which allows it to capture photons of selected 
wavelengths, while the reflectance of an object causes 
certain wavelengths to be re-radiated. A third character- 
istic affecting radiation is transmittance which is that 
property or characteristic of an object which allows the 
photons to pass unaltered through the object. 

Sources of Energy 

Sources of energy radiation may be generally cate- 
gorized into two types, broadband radiators and those 
producing line spectra or monochromatic radiation. For 
the most part, sensing in the visible and near-visible 
portions of the spectrum relies on broadband illumination 
or radiation. Solar radiation, a primary visible wave- 
length source, is broadband. It is only at the level of an 
individual molecule that line spectra radiation and absorp- 
tion become important. Laser illumination is of the single 
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wavelength variety and sensing a laser illuminated scene 
also requires special consideration. 

All objects at temperature above absolute zero (0°K) 
radiate energy. The theoretical limit of the spectral 
distribution of this radiation is described by Planck’s 
law ; 

_ 8irhv 2 

v 3 , hv/kT _ , 
c J (e - 1) 

where is the amount of radiant energy per unit frequency 
interval and v and T are the frequency and temperature, 
respectively. The constants c, k, and h are the speed of 
light, the Boltzman constant, and Planck’s constant. This 
radiation is characterized by an object which absorbs 
radiation equally as well as it radiates and with the 
same spectral weighting or distribution. Emission at the 
theoretical limit where a substance is a perfect absorber 
and emitter of energy is referred to as blackbody radiation. 

No real body is a perfect emitter, although some 
substances, such as lampblack or powdered graphite come 
close. They fail because they do not absorb or emit equal- 
ly over a wide spectral range. The ratio of the emittance 
of a real object, to that of a blackbody is called its 
emissivity and is a measure of the radiating efficiency of 
the object. The emission of real bodies behaves similarly 
to the blackbody laws in that the hotter the body, the more 
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the radiation and the shorter the dominant wavelength. 
However, the differences in the characteristic emission of 
a blackbody and that of a real body may be markedly differ- 
ent in certain spectral bands. 

Atmospheric Attenuation 

Since the nature of remote sensing is the measurement 
of objects at a distance, the effects of the intervening 
distance on the transmitted energy is important. In 
addition to the reduction of power density as an inverse 
function of the distance between the source and the detec- 
tor, the phenomena associated with atmospheric attenuation 
must be considered. Two major mechanisms are interacting 
with the energy flux density resulting in its attenuation. 
Specifically, these mechanisms are atmospheric absorption 
and atmospheric scattering. 

Beer's Law is a relationship which describes the 
change in flux density of a plane -parallel beam of mono- 
chromatic radiation as it penetrates a medium. This law 
is written: 

' y -ot.x 

h = h e x 

A 

where 1^ is the resultant flux density from the initial 
flux density I attenuated by the exponential factor with 

A 

parameters a^, the attenuation constant in nepers/meter, 
and x, the distance in meters through which the wave has 
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traveled. The distance x through which the wave must 
travel to be attenuated 37.3% or 1/e is referred to as the 
optical depth and has a value of 1 /a. . 

A major contribution to atmospheric attenuation, 
particularly in the visible and near-visible region of the 
electromagnetic spectrum, is molecular absorption by atmos- 
pheric constituents. It is important to realize that 
attenuation is related to both the rate of attenuation a, 

A 

and to the distance through which the energy must travel, 
x. In certain circumstances where the distance between 
sensor and target is small, the effects of the atmosphere 
will also be small. However, for satellite measurements 
where atmospheric absorption is high, the effects of the 
atmosphere can no longer be summarily neglected. 

In the ultraviolet portion of the spectrum, the 
nitrogen and oxygen absorption of radiation in the region 
from 0.13-0.22p is nearly complete. Radiation at 0.22- 
0.30u through an ozone layer is also completely blocked. 
Oxygen absorption bands occur at 1.06 and 1.27y in the 
infrared portion of the spectrum, while many absorption 
bands of water and carbon dioxide along with nitrous oxide, 
methane and ozone occur in the l-24y region. 

The atmosphere is completely opaque in many portions 
of the infrared spectrum. Absorption by atmospheric water 
vapor virtually closes the atmosphere from 25-lOOOp. Other 
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water vapor absorption characteristics of the atmosphere 
is found at 2.5mm, 5.0mm and 1.35cm. Despite these 
absorption characteristics of the atmosphere, many "win- 
dows" remain where the atmosphere is nearly transparent, 
and a significantly large percentage of energy is allowed 
to pass. Thus, the observation of the interaction of im- 
pending energy with a target is limited and in some sense 
is determined by the transmission characteristics of the 
atmosphere in the spectral band of interest. In Figure 
II- 2 is shown a graph of atmospheric transmission in 
various spectral bands. 

Aside from atmospheric absorption, there is also the 
problem of scattering by particles, water droplets and 
clouds within the atmosphere. Scattering is the process 
by which small particles suspended in a medium of different 
index of refraction diffuse a portion of the incident 
radiation. In most scattering (Raman being the major 
exception) no energy transformation results, only a change 
in the spatial distribution of the radiation results. 
Scattering is often a major cause of attenuation in the 
atmosphere . 

Scattering is a function of (1) the ratio of the 
particle diameter to the wavelength of the radiation; 

(2) the nature of the composition of the particles as 
characterized by their complex index of refraction; 



TRANSMISSION, % 
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(3) the particle size and number distribution; and (4) the 
shape of the scattering particles. 

When the size of the particle diameter is on the order 
of magnitude of about a tenth of a wavelength, scattering 
occurs in general accordance with the Rayleigh scattering 
law. According to the Rayleigh scattering law, scattering, 
hence attenuation, varies approximately with the fourth 
power of the wavelength. For larger particle sizes, on 
the order of magnitude of a wavelength, the scattering 
varies in a complex fashion as described by Mie theory. 

Mie scattering is characterized by re-radiation in the 
direction of the primary wave. Since the re-radiation of 
secondary waves is superimposed on the incident plane 
waves, not all of the energy is lost from transmission in 
the direction of the primary wave. A general relationship 
governing Mie scattering which is applicable to atmospheric 
aerosol in the visible and near infrared regions takes the 
form of S = C 2 A n , where the exponent n, which is a function 
of particle size, varies between -4 and a small positive 
number. At particle sizes ten times the wavelength, non- 
selective or wavelength independent scattering occurs. 

This particle size scattering effect generally serves to 
mark the somewhat diffuse upper bound of the realm of 
scattering and gives way to the laws of geometrical optics. 

Another form of scattering which is present, but 
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relatively small in comparison, is ’’Raman" scattering. In 
this type of scatter there is an inelastic collision 
between a photon and a molecule. The result is a change 
in the energy of the photons and consequently its wave- 
length. Thus the photon, in addition to being scattered, 
has a modified wavelength which is dependent upon the 
energy state of the molecule. The importance of the 
"Raman” affect is that it may be used to distinguish compo- 
sition and concentration of the molecular makeup of the 
atmosphere under special circumstances, but has a relative- 
ly miniscule affect when contrasted to other scatter 
mechanisms . 

Signatures 

Three basic energy interaction mechanisms are 
monitored by remote sensing instruments. These mechanisms 
include energy emission, absorption and reflection. 

Natural radiation is governed by Planck's law which 
describes the ideal emission from objects at temperatures 
above absolute zero (0°K). In Figure II-3 are shown the 
emission spectra for solar radiation at 3000°K and for 
emission from objects at 300°K, the average temperature of 
the earth. These curves show the dominance of solar 
radiation in the visible portion of the spectrum 
(0.4-0.75p) and the near infrared portion (0.75-3.0y). 






The radiation from the earth at 300°K becomes dominant in 
the thermal infrared portion of the spectrum (3. 0-14. Op). 

Solar radiation in the 0.4-3. Op portion of the spec- 
trum masks the radiation contribution of the earth even for 
low surface reflectances, therefore, it is this surface 
characteristic which most influences the energy measured. 

In this region, the reflectance and absorptance are for 
the most part complementary characteristics, and those 
substances or mixtures which have unique reflectance and 
absorptance spectra can be identified. The term ’’signature'’ 
has been commonly applied to spectral reflectance measure- 
ments of objects and characterizes a specific reflectance- 
absorptance relationship. In Figure II-4 are shown typical 
reflectance spectra for some natural objects. 

The radiation from earth targets at 300°K becomes 
dominant in the region of the spectrum from 3. 0-14. Op. It 
is in this region that the target temperature and emission 
characteristics strongly influence the measurements made. 
Measurements from the thermal infrared portion of the 
spectrum provide information about the temperature of an 
ob j ect . 

Beyond the thermal into the longer wavelengths lies 
the microwave region of the spectra. This region is more 
strongly influenced by the electromagnetic properties of 
volumes rather than strictly surface characteristics or 
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thermal perturbations, and may be sensed through passive 
radiometer measurements or by active radar techniques. 
Target signatures also exist in this spectral region which 
relate more strongly to the actual material composition of 
objects rather than to their surface characteristics alone. 
Although microwave sensing is rapidly developing, the 
remote sensing field has been largely dominated by visible 
and infrared sensing and these portions of the spectra 
have provided the bulk of the available measurements. 

Sensor Systems 

Sensor systems used in earth resources survey are 
varied in their techniques of energy measurement, the por- 
tion of the spectrum in which the energy is measured, and 
the ultimate format of the measurement presentation. 

Sensors may be described by their respective spectra, 
visible, infrared (reflective and thermal), and microwave; 
or they may be referenced by their format of presentation, 
either image or linear. They may be classed as either 
active or passive depending upon whether energy is radiated 
by the sensor and the return measured or whether the 
energy sources are external to the sensor and only their 
effects are measured. 

The most common sensor system is the camera. This 
imaging sensor is passive and photographically records 
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reflected energy from surfaces illuminated by an external 
source, typically the sun. Photographic sensing measures 
energy in the visible and reflective infrared portions of 
the spectrum (0.4-1. Op). This sensing is either broadband, 
where the energy density is integrated spectrally over the 
spectral range of the system (polychromatic black and 
white) , or the energy in spectral bands may be measured by 
selecting appropriate band pass optical filters and film 
combinations. Color film duplicates the spectral filter/ 
film combination wholly within the film. 

Other visible imaging systems belong to the multi- 
spectral scanner class of ins truments . Images are formed 
by multispectral scanner (MSS) systems by making measure- 
ments from an optical system which scans the area to be 
imaged. This instrument is passive and may be used to 
record the energy in several spectral bands simultaneously 
from a single target area. The multispectral scanner 
technique is also used to extend imaging into the reflec- 
tive and thermal infrared portions of the spectrum by an 
appropriate choice of system detectors. 

Sensor systems also exist which provide image sensing 
in the microwave region. These include both the active 
radar sensors and the passive microwave radiometer systems 
which provide a unique view of the earth. 
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Photography 

Photography is an important part of remote sensing 
activities. Not only are direct photographic recordings of 
remotely located objects frequently employed, but 
photography may also be the medium in which other sensor 
information is displayed. The principal advantage of the 
photographic process is its ability to quickly and accu- 
rately record radiant energy, and to do so economically. 

In order to utilize the photographic system to its full 
value, an understanding of the properties and limitations 
of the photographic process is necessary. 

Photographic film consists of a thin layer of base 
material, one or more emulsions and a subcoat which causes 
the emulsion to adhere to the base. The photographic 
emulsion is a specially prepared gelatin containing the 
photosensitive materials in suspension. Most emulsions 
also contain dyes to affect the photochemical process. 

Color film utilizes the basic photographic emulsion in 
three layers to produce the color image. Since the prin- 
cipal intention of color film is to replicate what the eye 
sees, only three dyes and three emulsions are necessary. 

An additional dye and emulsion layer would not be useful 
for normal viewing. 

In ordinary color film, the reception of blue light 
activates the appropriate dye system to replicate the blue 
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appearance by the processed film. Similarly, the received 
green light activates the dye system for replicating green 
and likewise for the red. In a color translation system 
such as infrared film, the received color activates a 
different dye system. In the case of infrared film, the 
blue light is filtered, the green light is made to produce 
blue on the film, yellow and red produce green and the 
received near-visible infrared produces the red excitance 
from the film. Thus, infrared reflection from an object 
is recorded and presented as red in the developed infrared 
film. 

Multispectral Scanner Systems 

Electromechanical scanner systems are imaging devices 
in which the scene to be mapped is physically scanned by a 
telescope system and the energy collected is converted to 
an electrical signal by a detector or detector array [ 18 ]. 
This method of image production has provided the majority 
of multispectral imaging available to date. The technique 
is reliable, flexible and proven. A considerable amount of 
development is being undertaken to replace the electro- 
mechanical systems with other electronically scanned 
devices such as the return beam vidicon, photodetector 
arrays, and the image dissector tube, however, these systems 
have not been adequately developed or proven to be fully 
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operational. The eventual replacement of the electro- 
mechanical scanner systems will not, however, change the 
basic sensing characteristics or image presentation modes 
produced by this general class of instruments. 

Scanner systems produce images by sweeping the 
detector image from side to side across the scene to be 
imaged while the scanner vehicle moves forward. This 
raster method of image scanning is graphically displayed 
in Figure 1 1 - 5 . The multispectral scanner technique also 
filters the monitored radiation into spectral bands prior 
to detection. Detection of the radiant energy in each of 
the spectral bands produces multiple images of the scene, 
which can be analyzed according to the individual spectral 
distribution of individual picture elements. The scanner 
system is an improvement over standard aerial photography 
in that both the spectral resolution and the absolute 
radiance calibration are more exact. With the magnetic 
tape recording of the signals produced by the scanner 
detection system, the data are in a format that is more 
amenable to automatic processing either by analog or by 
digital techniques. An additional advantage of the 
scanners is the fact that detectors are available for 
imaging outside the normal visible and near infrared ranges 
of photography. 
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Microwave Sensors 

Other imaging systems widely in use are the side- 
looking radar systems. These radar systems produce images 
by transmitting microwave pulses transverse to the direc- 
tion of flight. The radar return amplitude is measured 
as a function of time and is traced across a photographic 
film by an amplitude modulated glow lamp or similar device. 
This technique produces an image similar to the scanner 
system, however, the energy pulse transmitted is controlled 
so that an absolute measurement of the energy return is 
possible. The scanner system relies generally on the 
unknown radiative characteristics of the sun or of the 
terrain imaged. Synthetic aperture antenna techniques and 
other enhancement procedures may be utilized to obtain fine 
resolution and high image quality from the side-looking 
radar system. 

Electronically scanned array antennas may also be 
utilized with sensitive, low noise, broadband microwave 
receivers to produce microwave images. These images are 
produced as a result of measurement of the emissive 
characteristics of the terrain in selected microwave bands. 
Each of these microwave sensors provides additional 
knowledge of the terrain in a portion of the spectrum 
which has been relatively ignored. As the microwave 
characteristics of natural terrain features continue to be 
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investigated, great promise is shown for obtaining data 
with as high an information content as is currently 
obtained in the visible and near-visible region. Promise 
is also high for obtaining a better view of subsurface 
phenomena which cannot be directly detected by other than 
microwave techniques. 

ERTS 


The Earth Resources Technology Satellite was launched 
by the National Aeronautics and Space Administration in 
July 1972 to make repetitive, synoptic spectral measure- 
ments of the earth. The satellite was placed in an orbit 
such that nearly every point on the globe was visible to 
the satellite during an eighteen day cycle. The satellite 
system includes a four-band multispectral scanner imaging 
system which provides spectral images of the earth in 
green (,5-.6p), red (.6-.7p), and two reflective (.7-. 8 
and .8-1. Ip) infrared channels. These images are relayed 
to ground stations and are processed for wide dissemination 
in both image and digital computer compatible tape formats. 
For the first time, good synoptic and repetitive data are 
economically available routinely to a broad segment of the 
country. With the life span of ERTS-1 running short, the 
prospects for the launch of another ERTS series satellite 
become increasingly more promising. Apparently, this kind 
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of data, which holds the key to synoptic regional manage- 
ment of resources, is becoming increasingly and consis- 
tently available. It is also apparent from the amount of 
information which comprises an ERTS scene (2 million acre- 
sized picture elements) that traditional image interpretive 
techniques are inadequate and digital computer processing 
is the only viable method for data analysis. 

Analysis Systems 

Extensive analysis techniques have been developed for 
handling various forms of imagery. The large initial 
investment necessary to begin analyzing the more interest- 
ing problems offered by real data, as opposed to mathe- 
matical abstractions, accounts for the domination of this 
area of research by a few relatively large and well funded 
institutions. The analysis efforts have been concentrated 
largely in two major areas: (1) the extraction and en- 

hancement of monochromatic images, and (2) the assignment 
of specific categories according to multispectral signa- 
tures. Image analysis and enhancement techniques have 
been pursued largely by the Image Processing Laboratory 
at the Jet Propulsion Laboratory [19] and also at the 
University of Southern California [20]. Computer classi- 
fication from multispectral features has been introduced 
primarily at Purdue University [21], University of 
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Michigan [22], and University of Kansas [23]. Combinations 
of these techniques have been implemented on the NASA/JSC 
Digital Image Analysis System (DIAS) [24] . Software 
development in these areas has also taken place at the 
University of Maryland [7] and the University of 
Pennsylvania [25]. 

Recently, considerable emphasis has been directed at 
applying the techniques of computer graphics to the general 
areas of urban analysis systems, natural resource informa- 
tion systems and systems that concentrate on environmental 
problems [26]. Again the field is, to a great extent, 
dominated by the large facility, the computer- oriented 
designer and the sample solution, as opposed to an 
operational system. 

Although great strides have been made in the acquisi- 
tion, storage and automatic analysis of spatial and 
temporal data, the field is dominated by large computing 
machines and organizations, and intricate software systems. 
Yet according to the Rouse and Schell model [17], the 
benefit of the overall system must be measured by its 
effective use. Phillips [26] has observed, 

The most formidable obstacle to more 
effective use of computers in solving 
urban and environmental problems is 
the lack of effective "interf aces” 
between the user and the programs he 
uses. The interface problem is espe- 
cially acute because the user is not 
computer knowledgeable ... there has 
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been almost no progress either in 
presenting data in a form useful 
for problem solving, or in provid- 
ing functions to help transform 
data into useful information. 

Management Actions 

The resource management of large regional areas 
utilizing aerial imagery is quickly coming of age. Until 
recently, the technical capability of complex analysis of 
images was vested in only the largest groups which were 
able to afford the necessary computational equipment. This 
was complicated by the scarcity of the required timely 
images. With the advent of ERTS satellites and low cost 
minicomputer systems, the initial cost has become more 
realistic. Low resolution digital image displays capable 
of color enhancement are also becoming available at reason- 
able rates. The hardware components required to support 
regional management of resources seem readily available at 
this time. It is only the software aspects, the distilla- 
tion of image analysis techniques into an operationally 
usable system that is not as yet fully developed. Great 
numbers of analysis software packages have been written to 
address many aspects of the image enhancement and analysis 
mode, but all of these have been written requiring the 
support of great computing centers and learned practition- 
ers of the art of computer processing of image data. 
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It is important to remember that the utility of the 
whole remote sensing analysis system can be judged solely 
by the results of a user who is not computer knowledgeable. 
Nor can the user be generally expected to have an in-depth 
working knowledge of the physical characteristics of the 
measured energy, the engineering characteristics of the 
sensor systems, or the mathematical subtleties of the 
analysis systems commonly utilized in the remote sensing 
field. His expertise must be considered to be confined to 
the activity which he is best capable of performing, the 
analysis and planning related to urban, natural resource 
and environmental problem-solving. 

Graphic and cartographic methods for analysis of 
spatial data are required by the resource manager. The 
planner is concerned with extrapolating land usage and its 
associated environmental impact. Analysis of spatial 
comparisons is accomplished by a conceptual overlay of a 
variety of thematic maps which depict regions of arbitrary 
shape where attributes are homogeneously distributed. 
Variables are retrieved and assigned relative weightings 
corresponding to their desirability in the environmental 
design. Combining the weighted variables for a specific 
location in a functional form provides a value of the 
desirability of locations in terms of the environmental 
design. Other information of interest also includes the 
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determination of the areal extent or number of certain 
environmental or resource events, and data presentation is 
most usable in cartographic or statistical display format. 
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CHAPTER III 

INFORMATION SYSTEM DEFINITION AND SCOPE 
Introduction 

The application of remote sensor data to socio- 
economic problems is constrained by the rate of information 
flow across the data analyst-user interface. This rate of 
flow is impeded by the lack of awareness and understanding 
of each of these groups with the others’ capabilities and 
needs. A link is missing between the survey and analysis 
techniques of the remote sensing community and users. • 

Most current analysis systems are not developed to a 
fully operational state. Software still requires an 
analyst who fully understands the intricacy of its tech- 
niques and of complicated execution procedures, and 
processing systems are either highly complex, general 
purpose digital computers with equally complex operating 
procedures, or simple devices which rely primarily on 
heuristic analyses. Computer analysis techniques are 
fully operational only when they are developed to the 
point where the analysis steps are specified methodically 
and catalogued with use and parameter explanations. 

Complex data analysis techniques are required in a form 
compatible with users’ equally sophisticated, but solution- 
directed techniques. This void has been noted by an 
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applications panel convened by NASA [27] : 

The users comprising this application 
panel, as a result of their discussions 
and interaction with the different panel 
experts in the sensor technology, con- 
cluded that their current applications 
capabilities for sensor data lag behind 
the sensor technology. This is to say 
that these representative users acknowl- 
edged that their particular communities 
must for the most part improve their 
capability (in both numbers and sophis- 
tication) to assimilate, process, inter- 
pret, and apply the relevant sensor data 
which is essential to their programs. 

It follows therefore, that users, for the 
near term, must develop their own appli- 
cations capability. This is necessary 
before they can specify incremental 
additions to existing systems or the 
development of entirely new systems 
based on the applications needs. This 
does not imply, as a corollary, that 
research and development in the sensors 
field can mark time until the users 
catch up, but rather that additional 
R§D efforts for the near term must be 
directed toward enhancing existing sys- 
tems and improving the present quality 
of data reduction and processing. 

Personnel concerned with the interpretation and action 
on resource items available through an information system 
are non- technical , from a computer analysis point of view. 
For the most part, these people are ecologists, biologists 
and managerial personnel with relatively little or no 
experience in computer operations, computer programming, 
or the technologies involved in remote sensing data 
analyses. It is these individuals to whom this system 
development research is being directed, for these 



management functions for the most part have not been in- 
cluded in the specification of earth resources analysis 
systems. 

System Design Criteria 

Certain general characteristics are desirable in a 
resource information management system, especially in light 
of the expected system users. A specific system require- 
ment is a simplistic user interface which is interactive 
and self-explanatory in nature. A suitable dialogue to 
clearly state and explain the system analysis options is 
important, as well as the availability of ample error 
diagnostics and recovery procedures. It would appear that 
a system which implements the analysis capabilities 
required for data interaction, classification and enhance- 
ment is desirable. For the most part, analysis procedures 
with specified design parameters are required. A method- 
ology for the storage retrieval and handling of data, 
which is efficient, flexible and versatile, for the 
handling of source information forms the basic foundation 
of a resource information management and analysis system. 

For the most part, the general method for storage, 
retrieval and handling of sensor data currently used is on 
a per pixel (picture element) basis. This is convenient 
from a technologist's point of view for it is in this 
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format that the sensor data are available and it is in 
this format that image display is most easily handled. 
Although the per pixel basis is convenient and effective 
for the handling and display of data, this method falls 
short in the handling of the information derived from the 
interpreted image, whether it be manually or automatically 
interpreted . 

Resource information is obtainable by the interpreta- 
tion of the spectral distribution of energy reflected from 
the individual elements which comprise the imaged scene. 
This information is spatial and temporal in nature and 
available from varied sources. A data format amenable to 
the handling and comparison of multiple data sets which 
distinguishes between the point, line and area aspects is 
required by an information handling system. Additional 
capabilities for data transformation and data coordinate 
scale and location changes are appropriate. Suitable 
output product displays which relate directly to a decision 
making function of a user agency or which might supplement 

the manager's intuitive feel about a given situation are 
also important. 
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Data Base Representation 

Information Types 

An operational resource information system translates 
data from imagery, surveys and administration records and 
relates the data to a common geographic base. In a 
resource information base, integrity between three basic 
information storage and retrieval formats is desirable. 
These formats include information stored as : 1) point 

information; 2) lineal information, and; 3) areal informa- 
tion. 

Point information is that information which relates 
to a specific coordinate location or similarly, that 
coordinate point which locates the information. Standard 
sensor reflectance data are of this variety. A measurement 
of the spectral reflectance from a single ground resolution 
element is associated with the coordinate of that ground 
point. Also, any interpretation of the spectral data is 
associated with the point coordinate. Another association 
of this type might be the location of a specific water 
tower, for example. Here again is the interpretation/ 
coordinate association for a single ground point. 

Lineal information is characterized by the association 
of a single entity or data value for a number of ground 
resolution elements connected in a lineal fashion. 
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Examples of this type of association would include roads 
and area boundaries. This representation allows the 
aggregation of point information of the same value which 
may be represented by a larger more precise entity. 

Areal information is characterized by the association 
of data in an areal manner. Where adjacent ground cells 
have the same information associated with them, they may 
be accumulated under a single label and associated with 
the information. An example of this type of association 
might be all the resolution elements of a single cropped 
field (e.g. a cotton field) or might relate to a particular 
political association such as an area zoned for commercial 
development . 

Association of information under labeled aggregations 
of the type described allows for handling of large numbers 
of point/ inf ormation associations within a single aggre- 
gated unit. Thus, it not only reduces the storage required 
for the data, but allows for a more flexible interpretation 
and transformation of large volumes of data. Information 
stored under these separate types of labels is also easier 
to edit, expand and delete. The major problem in this 
labeled representation method is how to store, retrieve 
and analyze data with such widely varying data formats, and 
the definition of the specific formats of each of these 
data types. 
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Raster Format and Display 

A system of small area units based on geographic 
coordinates provides a flexible means of collecting, stor- 
ing, performing computer enhancement and decision opera- 
tions, and display of collected data. In this basic 
raster format, values are associated with a coordinate 
which is implied by a position or location in an array or 
string. Imagery is either collected or digitized generally 
by some form of raster scan. Many image sensors in the 
current earth resources survey inventory provide digitized 
data in this format. Digitization of photographic data 
also is accomplished in this format. Computer manipulation 
of the image elements as pixels in a random access stored 
array is extremely convenient, and image display of the 
image data by printer graphics or on color television units 
is greatly facilitated. It would seem that the availabil- 
ity of an image as picture element values in a random 
access array is the most suitable machine format. 

Only through. image interpretation analysis (either 
interactive, manual or automatic) do the sensor measure- 
ments become resource information, that is, knowledge 
required by an interpreter, analyst or manager which can 
be utilized in making a decision. The important charac- 
teristic is the location of a feature on the ground and 
its association with other similar, dissimilar or 
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aggregated features. In the sequential format, it is 
difficult to isolate individual units without reconstruct- 
ing a large portion of the data set in a random access 
storage medium. .This sequential digitization and storage 
method for all practical purposes precludes scale and 
coordinate rotations of the data for congruencing with 
similar data sets from the area. The storage method also 
is cumbersome in that all the data must be represented by 
individual data cells of the smallest integral size 
desired . 

Sensor measurements in and of themselves are useful 
only in as much as they provide additional knowledge about 
the location, relationship, condition, existence, etc. of 
the resource categories of concern. The termination of 
the digital analysis chain with the production of a raster 
format product, whether image or classified image, is 
unfortunate. It constrains the resource manager to the 
role of image interpreter, since the terminal product of 
almost every image analysis system is either the raster 
scan array stored on magnetic tape, magnetic disc, printed 
on a photograph or displayed on a CRT. 
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Boundary Representation 

Upon interpretation, the significance of individual 
picture element values is supplanted by the category 
assigned to their aggregation. A flexible, efficient 
methodology is important in the handling of these element 
aggregations. It is also desirable that this methodology 
be sufficiently flexible to handle information available 
from other sources relating information to spatial 
aggregations. Such a methodology can be built around a 
data format which aggregates ground points by defining 
the boundary which encloses or links them. This format 
can handle point, lineal and areal types of data in an 
efficient and flexible manner. 

A suitable data handling methodology for an informa- 
tion management system utilizes the boundary format, 
associating category values with the explicit coordinates 
of the boundary which links them. Utilization of paramet- 
ric functions which describe the boundary can minimize the 
amount of data stored. For example, if a piecewise linear 
boundary is assumed, only the coordinates of the line 
segments are saved as the information. In most cases, 
piecewise linear fits are adequate for data description, 
although quadratic or other parameter functions can also 
be utilized. 

Edge-tracing algorithms can automatically convert 
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raster information such as clas s if icat ion- interpreted 
images into the boundary format. Interactive conversion 
of interpreted information into the boundary format pro- 
vides the simplest approach to the format conversion. 

Figure I I I - 1 shows an area of Brazos County as imaged by 
the Earth Resources Technology Satellite and displayed by 
means of a line printer representation. The line printer 
image representation is typical of the raster format which 
is generally the analysis product. On the image product 
are shown interpretations of resource features which might 
be important to the area. These interpretations are 
representative of the point, line, and area information 
aggregations which are important. Piecewise linear 
boundaries are also indicated. 

The boundary format greatly reduces the amount of 
storage required for information maintenance and provides 
a flexibility unparalled by raster storage. As significant 
as the data reduction is the ease with which coordinate 
transformation may be enacted. With explicit coordinates , 
the transformations of translation, rotation and scaling 
become simple mathematical mappings involving no actual 
change in the elemental location of stored data. This 
format also facilitates the congruencing of data sets 
unlike either in scale or orientation. This boundary 
format representation is flexible, efficient, data- reducing 
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and convenient. 

The techniques appropriate for the conversion and 
handling of resource information in the boundary format 
are not particularly novel, however, the presentation of an 
information storage and retrieval format which provides a 
convenient and efficient methodology for transferring, 
handling, comparing and analyzing interpreted image sensor 
derived information is a significant step forward in the 
application of remote sensor data to the problems of 
resource management. 

Transformations and Data Handling 
Coordinate Transformations 

In a raster storage mode, the typical format for 
image remote sensor data, values are associated with a 
coordinate which is implied by the position of the value 
in an array or string. Coordinate transformation requires 
the repositioning of the value element in the storage 
array. Three basic coordinate transformations are re- 
quired for the handling of data. These are translation, 
rotation and scaling, which are extremely important when 
comparing data sets which differ in the time acquired 
sensor resolution, or format. Scaling and coordinate 
translation are not mechanically difficult for raster 
stored data, however, rotation is an extremely cumbersome 
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and time-consuming transformation. This is particularly 
true because coordinates are implied by the position at 
which the data are stored in an array. 

The coordinate transformations are no longer 
cumbersome procedures for data stored in the described 
boundary format. In this format the spatial coordinates 
of the data values are explicit and the storage of the 
data values is not constrained to a fixed location. 
Coordinate transformation of data in a boundary format 
relies solely on the mapping of the coordinate values 
associated with the information. With the explicit 
coordinate values, the transformation reverts simply to 
the application of linear coordinate transformation 
equations. The incorporation of an interactive coordinate 
transformation capability in an information management 
system is particularly important in congruencing data sets 
from differing origins or in locating information on 
imagery by overlay of familiar features. 

Data Transformation 

Resource inventory and control decisions are based 
particularly on observing the juxtaposition of favorable 
or unfavorable aspects of the environment. Under the 
"right" set of conditions, for example, a wildlife 
habitat may be suitable. A change in one of those condi- 
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tions may degrade the habitat to a marginal or even hostile 
one for the wildlife species under consideration. The 
kinds of information which are used, for example in 
determining a suitable habitat, are available from a 
variety of sources and may be presented in a variety of 
ways. The resource information management system which is 
used to aid in the making of environmental decisions 
requires the capability for logically comparing and 
combining decision elements into a single environmental 
quality value. Suitable mathematical and logical opera- 
tions can be interactively defined and utilized by the 
system operator. These operations include logical 
operations, conditional operators, and mathematical 
operators . 

Format Transformation 

An important aspect of a resource information 
management system is the ability to convert data, from one 
format to another. The boundary coordinate format 
described previously is suitable for an external format 
and a raster format is appropriate for display and image 
operations by an information management and analysis 
system. Conversion from an external boundary format to an 
internal raster format requires algorithms to expand 
point, line and area information. Interactive procedures 



53 


provide the simplest method for reversion from raster 
format back into a boundary format with the system operator 
providing the coordinate points defining the boundary and 
the information values associated with data entry. 

Analysis Operations 

The literature [6] describing the analysis of image 
data is resplendent with glowing reports of algorithms for 
spectral classification and image enhancement. These 
algorithms are the ’’answer” to the application of multi- 
spectral image sensing and provide meaningful inputs to 
the functional block labeled ’’user". These descriptions 
are just as conspicuously void of any description of how 
the image- derived data products can be applied and of any 
description of how the image formatted interpretations are 
available to the ’’user", unless he reverts to reproducing 
all the enhancement and classification results to dis- 
play them on his own system. These descriptions, however, • 
provide a comprehensive treatment of the various image 
analysis techniques which have been utilized to advantage 
in applying remote sensor data, and provide an adequate 
basis upon which algorithms may be selected and implemented 
in a resource information management and analysis system. 
Two fundamental analysis techniques required by the system 
are classification, and image enhancement. 
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Computer-Aided Classification 

Computer-aided classification is an algorithmic 
procedure for the computation and comparison of similarity 
measures between sensor measurements and the assignment 
of similar measurements to identified classes. A class 
or elemental signature is that set of measurements, 
characteristics or features which have been obtained from 
that class, or elements which describe the best knowledge 
obtainable by the sensor or sensors used to measure the 
object. The term "s ignature" is used to imply a unique 
association between the measurements and class elements, 
although the uniqueness of the association is not always 
supported. The primary assumption upon which classifica- 
tion is based is the assumption that measurements of 
dissimilar objects are not alike. The algorithm attempts 
to measure the similarities for classification purpose. 

There are two general approaches to the classification 
problem. These are supervised classification and unsuper- 
vised classification. Supervised classification is 
characterized by the a priori assignment of data elements to 
classes and the calculation from these assignments of 
typical class signatures against which unknown elements 
may be compared. The actual training procedure depends 
upon the actual classification algorithm utilized. 

The most commonly used classification algorithm is 
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based upon a maximum likelihood criterion. In this 
procedure, the form of probability density functions on the 
classes involved is assumed or determined. The measure of 
similarity between a data vector and a class is the 
a posteriori probability of the data point given the class 
which is most likely to have produced the point. The 
assumption is generally made that the distributions are 
multivariate Gaussian and the class sample means and 
covariances are used as estimators of the parameters of 
that distribution. This is the classification algorithm 
which forms the basis for the LARSYS series of classifi- 
cation algorithms developed at the Laboratory for Applica- 
tions of Remote Sensing at Purdue University [281. 

Unsupervised classification or clustering relies on 
the aggregation of similar data points and the a posteriori 
assignment of classes to these aggregates. The procedure 
may be single pass where data points are examined only 
once, or multiple pass where data points are examined 
several times in an attempt to improve upon the previously 
determined aggregations. Similarity measures are calcu- 
lated between the data points and the aggregations or 
clusters and the point is assumed to belong to the most 
similar cluster if the similarity measure is below some 
predetermined threshold. For this case, the cluster 
parameters are improved by the new data point and the 
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procedure continues. If cluster dispersion exceeds 
another threshold, the clusters may be split into two 
new clusters. If the data point is determined to be 
dissimilar from all defined clusters, a new cluster is 
formed by the data point. Differences in clustering 
algorithms are primarily due to differences in defined 
measures of similarity and in the procedures by which new 
clusters are formed or old clusters split. An attempt is 
made in some clustering algorithms to link clusters at the 
end of the procedure so that cluster shapes and sizes are 
not completely dependent upon the similarity measures and 
the thresholds chosen. Some versions of the LARSYS 
procedures use a combination of the two techniques by 
clustering within training set collections to improve upon 
the points selected for training and to give better 
estimates of the actual class signature parameters. 

Image Processing 

Certain image processing functions are desirable in 
an analysis system. These processing operations are 
utilized to improve the interpretabili ty of the imagery. 
Processing operations may be classed according to their 
mode of application: neighborhood, non-neighborhood and 

transformation operations. Neighborhood operations 
transform a data point according to the values of the 
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data points which surround it, that is the transformation 
is coordinate location dependent. An example of this 
class of operation is local averaging where a point is 
replaced by the average of the neighborhood to which it 
belongs. Non-neighborhood functions transform the data 
irrespective of its coordinate value or location. These 
operations involve scale changes in the image function 
such as grey-scale reversal or density slicing. Image 
transforms such as the Hadamard and Fourier transforms, 
deal with the mathematical mapping of the image space into 
another space. 

Several imaging processing functions are desirable 
in an image analysis system. Of particular importance 
are smoothing, edge enhancement and grey-scale transfor- 
mation operations. Smoothing and edge - enhancement may be 
accomplished either by neighborhood operations or by 
transformation. The simplest and most straightforward 
approach is the use of neighborhood operations. 

Communications 


Output Products 

Information management and analysis system output 
products are to a great extent limited by the types of 
peripheral output equipment available. Computer systems 
commonly include printer output and a CRT output capability 
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is desirable. General categories of output then are 
usually constrained to printer tabulation and graphics, 
and CRT image presentations. Other output to sequential 
file systems is appropriate for created or modified data 
files. 


Operator Communication 


The importance of an interacting man in an information 

system is emphasized by Martin [29]. 

...man must become the prime focus of 
system design. The computer is there 
to serve him, to obtain information 
for him, and to help him do his job. 

The ease with which he communicates 
with it will determine the extent to 
which he uses it. Whether or not he 
uses it powerfully will depend upon 
the man-machine language available 
to him and how well he is able to 
understand it. To be effective, 
systems will have to be designed from 
the outside in. The terminal or cas- 
ual operator instead of being a 
peripheral consideration will become 
the tail that wags the whole dog. 

Design of the dialogue for operator communications 
requires consideration of the level of training of the 
expected operators of the system. Although these 
individuals are intelligent, their exposure to computer 
languages and processing operations is generally limited. 
The operators for the most part are casual operators as 
opposed to operators dedicated to system operation. These 
casual operators need to use an information and analysis 
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system at some irregular frequency and for the most part 
cannot be expected to remember complex mnemonic entry nor 
the full sequence of operations or capabilities of the 
system. 

In view of these operator considerations, the 
computer-initiated dialogue is preferred, refreshing the 
operator on the full capability of the system whenever 
required. Certain functions of the system require entry 
of programming- like statements, in particular for coor- 
dinate transformation or for data comparison and decision, 
however, this interaction should be kept simple and to a 
minimum. 

Another procedure for communication with the operator 
described is menu selection. Fixed responses are expected 
in particular for specification of the processing tech- 
niques. Display of the various options in a menu format 
with a method for selection of the desired technique offers 
a reasonable approach to operator communication. Flexible 
but simple and informative operator communications are 
necessary for an effective interactive system. 

Error Diagnostics 

Terminal operators usually make errors. It is essen- 
tial to protect the computer system and its data from 
these errors, as far as is possible, and to do this 
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without antagonizing the operators. A factor that is 
strongly in favor of real-time systems is that in an appro- 
priately designed dialogue, most operator errors can be 
caught in real-time. The mistake or discrepancy is 
detected and corrected on the spot. The effectiveness of 
error detection is highly dependent on the design of the 
man-machine dialogue. 

Psychological considerations are important in the 
design of a man-machine dialogue. It is desirable to 
steer a course between boredom and bewilderment. Studies 
of the human learning process have shown that positive 
response to correct actions and admonition of incorrect 
actions within seconds give by far the best reinforcement 
[30]. The conscientious operator learns from the response 
to his incorrect actions. The response, however, need 
not be overly abrupt: a split second error response amid 

thought is jarring and rude. The error response suitably 
is delayed until the operator is permitted to finish his 
thought. One danger in giving real-time error responses 
is a tendency for an operator to become careless expecting 
the machine to catch all of his mistakes. 

The requirement for error detection and diagnostic 
response is apparent. However, a dialogue critically 
examined in light of both motivational and technical con- 
siderations is more effective. 
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IRIMAS--An Interactive Resource Information Management 

and Analysis System 

The specifications outlined above are embodied in an 
interactive resource information and analysis system-- 
IRIMAS developed on a Texas Instruments 980A minicomputer 
system. This is an interactive software system configured 
for simplistic analysis of resource information derived 
from remote sensor data interpretation. It is also suited 
for inclusion of resource information from other informa- 
tion sources. The salient features of the system include 
a data handling system for the storage, retrieval and 
analysis of resource information in the boundary format 
described. The data are encoded for extended -storage in 
data strings where the boundaries are defined by piecewise 
linear curve fits to the coordinates. System flexibility 
is maintained through the utilization of syntax definitions 
and syntax interpretation algorithms. 

Syntax algorithms are utilized for: (1) conversion 

of the external information character string into an 
internal format suitable for further analysis; (2) the 
interactive definition of data- and coordinate- trans forma- 
tions in programming- like statements, and; (3) the inter- 
active operation interface to the information system for 
the specification of system analysis functions. 

Algorithms for the analysis of remote sensing and 
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information data are implemented to include the following 
capabilities: (1) CRT and printer map display, (2) maxi- 

mum likelihood classification, (3) image enhancement with 
smoothing, edge enhancement, and single pass clustering, 
and (4) computation of data statistics. 

An interactive operator interface is implemented to 
allow simplistic operation of IRIMAS by specification of 
function procedures for standardized data analysis. 

In this chapter are discussed the general requirements 
for a computer system to perform interactive analysis of 
resource information and image interpretation. A mini- 
computer system is developed using these general guidelines 
and specifications to address the problem of interfacing 
image data available from remote sensing sources to the 
resource management function. The specific components 
and implementation of this computer system, IRIMAS are 
subsequently described in following chapters. 



63 


CHAPTER IV 

IRIMAS DATA HANDLING AND TRANSFORMATION COMPONENTS 

The general design goals for an interactive informa- 
tion and analysis system which aids information management 
and image interpretation by a resource planning function 
are described in previous chapters. An Interactive 
Resource Information Management and Analysis System, 

IRIMAS, developed with these design goals is implemented 
on a Texas Instruments 980A minicomputer system. Although 
the implementation is restricted to use on the TI 980A, 
the techniques employed are applicable to any general 
purpose computer system with suitable display peripherals. 
Subsequent discussion describes the techniques applied in 
IRIMAS and the specific implementation employed. 

The more significant features of the general informa- 
tion management system described and of IRIMAS in particu- 
lar are concerned with the encoding, decoding and handling 
of point, line and area information. In IRIMAS, the 
boundary format previously discussed is used for informa- 
tion storage and retrieval. The flexibility required by 
the system is met through the use of syntactic definitions 
and algorithms for syntax interpretation. Configured in 
this way, modification of the interactive input and output 
procedures and changes in the specific format utilized for 
information representation mostly involve syntax definition 
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change. These changes are localized to specific sub- 
routines which are easily modified; the major logic and 
processing components of the system remain unchanged. 

In this chapter are discussed the background required 
for the understanding of formal language procedures and 
the software devices required for interpretation. The 
specific components of IRIMAS required for interactive 
transformation of data and coordinates are discussed, as 
well as those components required for interpretation of 
the information storage format. System algorithms and 
software components used in expanding the boundary 
coordinate information representation into the raster or 
grid data arrays is also covered. 

Formal Languages 

Formal language theory has developed an abstract 
terminology for the modeling and representation of 
information structures which can be represented by a set 
of primitive elements together with a set of rules for 
combining those primitives. The formalization of languages 
and their representation has allowed the development of 
general procedures for the recognition of sentences in a 
suitably defined language without undue concern over a 
specific language definition. A subclass of all languages 
is the set of languages defined by phrase- s tructured 
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grammars which have been widely used in the definition of 
formal computer programming languages such as FORTRAN, 
ALGOL, etc. The translation of these programming 
languages is done mechanically by compiler software which 
implements the general syntax-driven parsing procedures. 
These language techniques have been used to define a data 
structure for the data maintained by the IRIMAS information 
system, and for the definition of operator language and 
command, and interpretation. 

Language Terminology 

The terminology common to formal language specifica- 
tion is defined in references [ 31] and [32] . An a lphabet 
is any finite set of symbols used to converse in a lan- 
guage. In the English language, this would correspond to 
the well-known 26 letter alphabet plus those symbols or 
punctuation which might also be appropriate. A sentence 
is any string (or sequence) of finite length, composed of 
symbols from the alphabet. A language defines a set of 
sentences over an alphabet domain. In the formalization 
of a language, the need to recognize what alphabet strings 
belong to the language, where the number of sentences in a 
language is infinite, becomes a problem. For certain 
languages, the sentences of the language can be generated 
by a procedure which is formalized by the concept of a 
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grammar. A grammar is the formal definition of those 
sentences which generate a language. The structural 
definition of the language in terms of a grammar is called 
a syntax of the language. 

These terms and their relationship might best be 
demonstrated by a sentence from the English language. As 
an example take the sentence, "The black cat ran away." 

In terms of English grammar, this sentence is composed of 
a noun phrase, "the black cat", and a verb phrase "ran 
away". The noun phrase consists of the modifiers "the" 
and "black", and the noun "cat", while the verb phrase 
consists of the verb "ran" and the adverb "away". The 
grammar for this sentence is defined in terms of (1) syn- 
tactic categories from which strings of words could be 
derived; (2) variables or non- terminals ; (3) terminals or 
words; and (4) productions, which are the relations that 
exist between various strings of variables and terminals. 

A grammar also specifies a non-terminal "sentence" or 
"start" symbol which generates exactly those strings or 
terminals that are in a language. 

The presentation of the example sentence in terms of 
a formal grammar and its notation is as follows: 

<sentence> : : = <noun phrasexverb phrase> 

<noun phrase> : : =<modif ierxnoun phrase> | <noun> 

<modif ier> : : = The | Black 
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<noun> : : = Cat 

<verb phrase> : : =< verb><adverb> | <verb> 

<verb> : : =<ran> 

<adverb> : : =<away> 

The grammar consists of non- terminals or variables 
such as <noun phrase>, <noun>, <adjective>; terminals or 
words such as "the," "cat," "black," and productions are 
those relations that exist between various strings of 
variables and terminals. The symbols stands for 

"is defined to be", "|" stands for "or" and the angular 
brackets are used to enclose each name of a phrase or 
variable. The formal grammar also requires a starting 
symbol which in this case is the variable <sentence>. 

An additional method for representing a grammar is 

in terms of a syntax table. A syntax table is the tabular 

\ 

form extensively utilized by the information syntax 
algorithms. A syntax table for the above sentence may be 
constructed as shown in Table IV-1. 

Syntax Driver 

A general purpose driver may be formulated to step 
through syntax tables to determine if a sentence belongs 
to the defined grammar. The driver will try to match a 
sentence like "The black cat ran away" to the syntax 
tables. In doing this, it must keep track of the current 
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Example Syntax Table 


s entence 

equation 
equation 
end of table 

noun phrase 
verb phrase 

noun phrase 

equation 

or 

equation 
equation 
end of table 

noun 

adj ective 
noun phrase 

adj ective 

literal 

or 

literal 
end of table 

The 

Black 

noun 

literal 
end of table 

Cat 

verb phrase 

equation 
equation 
end of table 

verb 

adverb 

verb 

literal 
end of table 

Ran 

adverb 

literal 
end of table 

Away 
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syntax position. 

The sentence is held, internal to the computer, as a 
string of words. An input cursor is used by the driver 
to position itself in the sentence. When processing 
begins, the cursor will point to the first word in the 
sentence and as processing continues, the input cursor 
will advance. 

As the sentence is processed, sentence words are 
matched against items in the syntax tables. The driver 
will keep track of its current syntax table position as it 
sequences through the table. With each syntax table 
position is associated a syntax position type and a posi- 
tion value. The syntax position type designation indicates 
whether the position value is a variable (equation), a 
terminal value (literal),- or the type dictates a table 
delimiting function (or, end of table). Two additional 
types of syntax position are also utilized. The position 
type "repeat” indicates that at this point in the syntax 
table the non-terminal value is considered to be repeated 
zero to an arbitrary number of times. Syntax position 
type "subroutine" indicates that the position value is to 
be output to the interpreting subroutine as an operation 
code. The position of an output string maintained by the 
interpreting subroutine is also available to the driver. 

The driver maintains push-down stack storage of the 
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input position, the syntax position and the output posi- 
tion. It keeps a series of these pointers in the stack 
as it works its way through the table for each equation. 
Then, when the driver is ready to return through the 
equations, it has a memory of its path. Having taken a 
wrong path through the syntax equations, the driver also 
has a memory of input and output cursors and thus can re- 
store the input and output cursors from the stack to try 
an alternate path. 

In Figure IV- 1 is shown a flow diagram of a general 
purpose syntax interpreting algorithm. Syntax interpre- 
tation is handled by four components of the syntax driver 
corresponding to the syntax table position types: 

(1) equation/repeat , (2) literal, (3) or/end of table, 
and (4) subroutine call. As the type of current syntax 
position is determined the appropriate software component 
is utilized to further process the input string. The 
syntax driver outputs operation codes by way of the sub- 
routine call and halts with a "valid syntax" decision if 
the input string exactly matches the syntax table, or an 
"invalid syntax" decision if the two do not match. For a 
detailed description of the syntax driven algorithm see 
Gauthier and Ponto [34]. 
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Figure IV- 1. Logical Flow Diagram for 
Syntax Interpretation 






72 



Figure IV-1. Logical Flow Diagram for 

Syntax Interpretation (cont.) 
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Figure IV- 1. Logical Flow Diagram for 

Syntax Interpretation (cont.) 




XRIMAS Information Syntax 


A flexible resource information data system must 
handle data files containing point, lineal and areal infor- 
mation entries. These entries must be referenced to a 
common geographical coordinate base and must be flexible 
enough in their definition to include a broad variety of 
numerical, literal and logical information. A data base 
with such varied and flexible entries can best be built 
and interpreted if syntactic representation and parsing 
algorithms are used. This notational representation 
facilitates the definition of a hierarchial structure and 
allows mechanical interpretation or parsing of an arbitrar- 
ily defined language within a specific class of grammars. 

The format of the IRIMAS data base is defined in 
terms of character strings which are collected into phrases 
and sentences similar to a spoken language. Each sentence 
has a different functional purpose in the precise defini- 
tion of the information type and its value. Sentences in 
this data base definition have the following purposes: 

(1) data set identification; (2) attribute definition; and 
(3) data definition. Data definition statements provide 
point, lineal and areal data entry. Additionally, a grid 
or raster format has also been included as a part of the 
information system format. 

The data set is of a sequential nature and the first 
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item is a data set identification and description label. 
Items included in this label are the file name, geographic 
coordinate references and a narrative description of the 
data file to include such items as origin and previous 
processing, together with other descriptive material. 
Optional data set inclusions following the description 
label are "sentences" describing and defining the list of 
attributes associated with the data file. Numerical codes 
used to reference lengthy data cell descriptions would be 
defined in these data set entries. The third type of data 
set entries are the sentences defining the data. These 
entries would include a data entry label which describes 
the entry name and the type of entry (point, lineal, areal, 
or grid) , a coordinate definition referenced to the 
coordinate definition included as part of the data set 
label, and the values of attributes applying to this 
specific data entry. These attribute values could either 
be quantitative or descriptive. Other items included in 
the data set would be delimiting punctuation between 
sentence entry items. 

Block Label 

The first entry in a block of data entries is a 
Block Label which contains identification and descriptive 
information about the entries to follow. This specific 
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entry is free form, where specific elements are position- 
ally located. That is, items in the entry must follow a 

fixed sequence, but their lengths and absolute location in 
the data record are unconstrained syntactically. All 
characters in the entry are alphanumeric. An entry 
identifier, the first byte of the entry, is arbitrarily 
assigned the ASCII value of Thus, the first item in 

the sentence is the sentence identifier. Subsequent items 
are separated by either spaces or commas. The second item 
to be included is the data set name. This is the label by 
which the block is uniquely identified and logs may be 
kept relating to the origin of the information contained 
in the data set. The third and fourth entry items are 
the geographical coordinate reference for the block. 

Their exact natures are undefined here but longitude- 
latitude might be a typical entry item. 

A narrative description of the data set may follow 
the coordinate values. This description would include any 
information pertinent to the data set which would be 
necessary at analysis time. The last item in the entry is 
a delimiter code arbitrarily selected to be an ASCII "/". 
This delimiter is used to separate sentence entries into 
the data set much like a period in a regular language. 
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Attribute Definition 

The second type of entry, the Attribute Definition , 
allowed in the data set provides for the more precise 
definition of data attributes. This entry has no opera- 
tional significance and is purely descriptive. Its 
purpose is to provide a more complete description of the 
values assigned to an attribute variable. Items within 
this entry are delimited by colons (:) and commas. Spaces 
are recognized only as valid characters. The first item 
in this entry, like all other entries, is a defining 
operation code, selected in this case to be an ASCII 
The second item is the attribute name to be described. 

This name is separated from the subsequent descriptive 
information by a colon. The additional items in the entry 
are used to describe values assigned to the attribute in 
later data entries (e.g. 0 = no classification selected, 

1 = cotton). These descriptive items are delimited by 
commas and the attribute definition entry is terminated by 
the delimiter 

Data Entries 

The third type of entry allowed in the data set is 
the Data Definition statement. This sentence consists of 
several phrases each identified by a characteristic opera- 
tion code and each phrase includes several items. A data 
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label phrase is used to initiate a data definition entry. 
This phrase includes an ASCII operation code as the 
first item, which designates a standard type data entry. 
The next item is an alphanumeric "P M , M L M , "A", or "G” 
to define the data type. The letter "P" designates a 
point entry, "L" a lineal entry, "A" an areal entry and 
"G" a grid entry. The final item in the data label is 
the unique name associated with the entry. No delimiters 
are recognized in this phrase, which is terminated by the 
following phrase operation code. Once the data definition 
entry is opened by the data label phrase, the order of 
additional entry phrases is not restricted. 

A second type of phrase for the data entry is the 
coordinate definition. This phrase begins with the opera- 
tion code ASCII M % M . A single occurrence of this phrase 
is sufficient for a point data entry; multiple coordinate 
entries, however, are required for both lineal and areal 
entries. Lineal entries are defined by the sequence of 
coordinates derived by approximating the curve with 
straight line segments. The coordinate entries are the 
endpoints on the curve of those lines. Areal entry 
coordinates are similarly derived only the curve closes 
on itself. The area defined is as that enclosed by the 
curve and the attributes entries refer to that area. 
Attributes of a lineal entry describe the curve. 
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The second entry in the phrase is the coordinate 
variable (either "X” or "Y") followed by the coordinate 
value. These items are followed by the second variable 
name and the coordinate value where applicable. In a 
differential mode, these additional items are optional. 

The omission of a coordinate variable definition in the 
differential mode is taken to mean no change in that 
variable value from the proceding coordinate entry. 

The mode of the coordinate insertion is the fourth 
item in the phrase. An alphameric "A M designates 
absolute and an alphameric "D" designates the differential 
mode. Note that each data entry must have at least one 
absolute coordinate definition as the first coordinate 
definition. 

Two types of data value- defining phrases are 
specified. The first beginning with an operation code 
n e n includes attribute values which are descriptive in 
nature. The second item in the phrase identifies the 
coordinate by its position in the block definition 
statement. Alphamerics are used for this identification. 
The letter "A” designates the first attribute "B" the 
second attribute, etc. This attribute designation is 
followed by the descriptive attribute value. For the 
quantitative attribute specification, the phrase begins 
with the operation code ,f # M followed by the attribute 
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designator as above. The next item in the phrase is the 
numeric string defining the value of the attribute. These 
phrases are terminated either by the beginning of another 
phrase within the entry or the entry delimiter 

The IRIMAS data entry statement also provides for 
gridded format data entry. Entry statements for gridded 
data utilize the coordinate entry format to specify the 
upper left-hand corner and the lower right-hand corner of 
the rectangular data block entered. The first attribute 
"A" is the number of 8-bit bytes in each picture element 
in the gridded array. After the data entry statement, the 
data are entered in a sequential manner with x-coordinate 
values changing most rapidly. The data entry statement 
terminates after all the specified data has been input. 

Formal Syntax Definition 

The statement definitions described above may be 
more formally represented in terms of syntactic defini- 
tions : 

<entry> : :=<block label> | <attribute> | <data> 

<block label>::= £<name><block coord. ><description> 
<attribute> : : = § <attribute name>/<description> 
<data>: :=<data label><coordinate><data> 

<name>: : =/<alphanumerics> 

<block coordinates> : : =<longitude><latitude> 
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<longitude> : : =< degrees>^<minutes>j_<seconds> 

<latitude> : : = <degrees>j_<minutes>j_<seconds> 

<degrees> : : =/<numeric> 

<minutes>: : =/<numeric> 

<seconds> : : =/<numeric> 

<alphanumerics> : : =/<alphamerics>/<numerics> 
<description> : : =/< alphanumerics> 

<data label> : : =*< type><name> 

<type> : : - A | P | L \ G 

<coordinates> : := %<mode>/ [<variable>/<numerics> ] 
<variables> : : = X|Y 

<data> : :-< descrip tive> | <quantitative> 

<descriptive> : := @ <attribute><alphanumerics> 
<attribute>: : =<alphamerics> 

<quantitative> : := #<attribute>/<numerics > 

<alphameric> : : = A | B |£| D. . . | 

<numerics> : : = + | • • |£ 

Machine Format 

The data format specified for information storage and 
retrieval is inconvenient for manipulation and analysis by 
computer software. Consequently, syntax driver parsing 
algorithms are used by IRIMAS to convert the system data 
language into a convenient format for machine analysis. 
These parsing algorithms convert the syntax data into the 
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constituent elements of alphanumeric descriptions, 
coordinates, data type, attribute values, and also set 
attribute flags to indicate that specific attributes are 
related to descriptive text rather than quantitative value. 

The data for each entry element are maintained in 
several arrays. The array containing the alphanumeric 
descriptors contains the entry element name as the first 
string in the array. Additional strings in the array are 
produced by those attribute entries which have descriptive 
values rather than quantitative values. The length of the 
string is stored as the first word-sized element in the 
string and the string is word- aligned. Coordinate data is 
stored as x-coordinate and y-coordinate pairs in a stack 
array. The first element in the array is the number of 
entries in the stack. The coordinate pairs are stored as 
absolute values referenced to the data set beginning 
coordinates rather than in the differential mode which is 
also permissible in the data entry format. 

The attribute array contains several kinds of infor- 
mation. Attribute numbers are fixed and each attribute 
entry is either an attribute value for those attributes 
with quantitative values, or the position contains a 
pointer to the string in the alphanumeric descriptor array 
which corresponds to a description value of the attribute. 
In order that there can be no confusion as to whether an 
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entry is a value or a pointer, a bit is set to one in the 
first two words of the actual attribute array. The bits 
of the double word correspond to the attribute position in 
the array with the high order bit corresponding to the 
first specified attribute. The third coordinate in the 
array indicates the type of data entry with point, linear, 
area and grid entries corresponding to type values of 0-3, 
respectively. This iri essence leaves the fourth position 
attribute, designated '’A' 1 in the data entry format. 

Grid format entries from the data set utilize the 
coordinate array to store the coordinate extents of the 
grid input, and utilize the first attribute value to 
specify the number of bytes corresponding to each coordi- 
nate point entry. Specific quantitative values input in 
the grid format are maintained in an additional buffer 
area reserved for grid data. It is in this machine format 
that the data entered into the IRIMAS software system is 
further analyzed. 

Data Syntax Conversion Subroutine 

The IRIMAS system uses a conversion subroutine 
SUBR.l, together with the syntax driver algorithm to 
convert the data as interpreted by the driver into the 
machine data format. Conversion is controlled through a 
set of eight operation codes to the subroutine which are 
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embedded in the syntax tables defining the information 
system data format (see Table IV-2). Additionally, the 
current position of the input cursor is provided to the 
syntax driver to allow correction of data for improper 
traces through the syntax table. 

Several buffer areas are maintained within the con- 
version subroutine software. Buffers for current x and y 
coordinate values (CV.X, CV.Y) are maintained as well as 
a buffer location (TEST) which is used to hold characters 
from the input string for comparison. A pointer to the 
output buffer containing data alphanumeric descriptions is 
maintained which points to the beginning of the character 
string associated with the current syntax variable being 
processed by the conversion subroutine. The output cursor 
generally points to the location of the last character in 
the output buffer. An additional buffer area (VALUE) is 
provided which holds the address pointer to the location 
of the current attribute or coordinate value being 
processed. 

Conversion subroutine operation codes, their syntax 
implications and subroutine processing are described in 
Table IV- 3. In Figure IV-2 is shown the logical flow 
chart for the conversion subroutine. 

When the syntax-driven algorithm is utilized together 
with the appropriate syntax table and the conversion 



TABLE IV- 2 


Information Syntax 


Entry 

Subroutine 

Repeat 

Equation 

Repeat 

End of Table 

Cr. Atr. 

Equation 

Or 

Equation 
End of Table 

Attribute 

Equation 

Or 

Equation 
End of Table 

Label 

Literal 

Equation 

Subroutine 

Equation 

Subroutine 

End of Table 

Type 

Subroutine 

Literal 

Or 

Literal 

Or 

Literal 

Or 

Literal 
End of Table 

Name 

Repeat 

End of Table 

Alphanumeric 

Equation 

Or 

Equation 
End of Table 


0 

Delimit 
Label 
Cr. Atr. 


Coordinate 

Attribute 


Data- Alpha 
Data-Numeric 


* 

Type 

1 

Name 

Alpha 

3 

A 

P 

L 

G 

Alphanumeric 

Alpha 

Numeral 
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Alpha 


Numeric 


Coordinate 


Coord Value 


Mode 


TABLE IV- 2 (cont.) 


Subroutine 

Literal 

Or 

Literal 

Or’ 

Literal 

Or 

Literal 

Or 

Literal 

Or 

Literal 
End of Table 

Subroutine 

Literal 

Or 

Literal 

Or 

Literal 

Or 

Literal 

Or* 

Literal 
End of Table 

Literal 
Repeat 
Equation 
Subroutine 
End of Table 

Equation 
Subroutine 
Repeat 
Subroutine 
End of Table 

Subroutine 

Literal 

Or 

Literal 
End of Table 


3 

A 

B 


Z 


IT If 


3 


1 


9 


% 

Coord value 
Mode 
6 


Value 

4 

Numeric 

5 


3 

A 

D 



TABLE IV- 2 (cont.) 


Variable 


Data- Alpha 


Data- Numeric 


Subroutine 

3 

Literal 

Or 

X 

Literal 
End of Table 

Y 

Literal 

@ 

Equation 

Alpha 

Subroutine 

7 

Repeat 

Alphanumeric 

Subroutine 
End of Table 

2 

Literal 

# 

Equation 

Alpha 

Subroutine 

7 

Repeat 

Numeric 

Subroutine 
End of Table 

5 

Literal 
End of Table 

/ 


Delimit 



Table IV- 3 


SUBR.l Operation Codes 


Operation Code 0 

This operation allows the subroutine to 
initialize the machine data arrays and 
conversion subroutine parameters. 

Operation Code 1 

This code indicates that the input cursor 
points to the alphameric code for the data 
entry type. The conversion subroutine sets 
the type value in the machine buffers with 
values 0-3 assigned to point, line, area 
and grid data, respectively. 

Operation Code 2 

An alphanumeric string is terminated. 
Initially, V-ALttE points to itself which 
indicates to the subroutine that the first 
alphanumeric string is the name of the data 
entry. If this is the case, the subroutine 
returns to the driver. If VALUE does not 
point to itself, the alphanumeric string 
processed represents the descriptive value 
of an attribute in the attribute array. 

In this case, the address of the first 
character in the string which has been 
moved to the output buffer is placed in 
the machine location assigned to the 
attribute. For alphanumeric description 
valued attributes, the attribute location 
in the attribute array contains an address 
pointer to the description string. 

Operation Code 3 

The literal at the input cursor is moved 
to the output buffer and the output cursor 
is incremented. 

Operation Code 4 

The syntax driver has identified a coor- 
dinate variable at the input cursor. The 
pointer VALUE is set to the address of the 
coordinate buffer corresponding to the 
literal at the input cursor, "X" for 
x- coordinates and "Y" for y- coordinates . 

The pointer POINT is set to the current 



Table IV- 3 continued 


value of the output cursor. 

Operation Code 5 

This operation code indicates that the 
character string moved to the output 
buffer is an ASCII coded decimal value. 
The subroutine converts this decimal 
integer into a binary integer and 
stores the binary number at the address 
indicated by the pointer VALUE. 

Operation Code 6 

The end of a coordinate data entry has 
been detected and the input cursor 
points to the coded value for the mode 
of the coordinate data point. If the 
coordinate mode is absolute, the values 
in the coordinate buffers are stored in 
the coordinate values array. If the 
coordinate mode is differential, the 
previous coordinates are added to the 
values in the coordinate buffers and 
the resultant coordinate values are 
stored in the coordinate stack. 

Operation Code 7 

This operation code indicated that an 
attribute variable has been indicated 
and the attribute code is at the input 
corner. The subroutine sets the VALUE 
pointer to the address of the specific 
attribute in the attribute array and 
sets the pointer POINT to the value of 
the output cursor. 
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Figure IV- 2 , Logical Flow Diagram for SUBR. 1 
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Figure IV- 2. Logical Flow Diagram for SUBR.l (cont.) 
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Figure IV- 2. Logical Flow Diagram for SUBR.l (cont.) 
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subroutine, input character text following the syntax 
definitions the equivalent machine variables are set 
according to the information input. Processing of the data 
in machine format may proceed from this point. 

Data Transformation 


Definition 

A desirable feature incorporated into IRIMAS is the 
capability of combining weighted variables together in a 
functional form and the comparing of the variables and 
derived values. In order that these data transformations 
are programmable real-time, an interactive programming- like 
operator interface has been established. Program- like 
statements are input by a system operator and are inter- 
preted by a syntax driver, syntax table and conversion 
subroutine. The output string from the conversion sub- 
routine is interpreted to perform simple mathematical and 
logical functions on the machine format data. 

Three specific sets of variables may be referenced: 

(1) a set of dummy variables; (2) the output variables 
which are kept separate from the input variables; and 
(3) the attribute values themselves. Each variable within 
the range of these arrays may be utilized by the inter- 
pretation scheme. The following types of data transfor- 
mation, and logical operators are allowed: 
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Logical 
. EQ. 
.NE. 
.GT. 

. LT. 

. GE. 

. LE . 


Arithmetic 
- subtraction 
+ addition 
* multiplication 
/ division 
** Exponentation 


Arithmetic subexpressions are also permissible in the 
data system as well as conditional statements like 


IF (A1 = 0.0) THEN (A2 = 3.0). 

A limited number of variables are allowed for this 
interpretative phase. These variables are referenced by 
a prefix letter and a suffix number. Attribute variables 
are designated with the prefix "A" and twenty seven may be 
referenced (A0-A26). The first attribute variable (AO) 
designates the data entry type (point, area. .. etc. ) . A 
set of dummy variables are available for programming 
flexibility. These are designated with a prefix "D" and 
range from D0-D4. Similarly, output variables are 
designated with a previx "0" and range from 00-04. 
Statements allowed are FORTRAN- like equations, however, no 
intrinsic functions (sine, cosine, etc.) are supported. 
Neither are any branch or input/output operations permit- 
ted. An example of a permissible set of statements is as 
follows : 
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00 - 0 

DO * 3. 0*A1 + 0.9 * A2 + 1*A3 
D1 - A24 - A15/2 
IF (Dl.GT.DO) THEN (00 = 255) 

END.. 

The programming statements are defined by a syntax 
table and are converted into Polish algebra strings which 
are interpreted to effect the mathematical operations. 
Polish notation is described in [38], Essentially, the 
technique presents numerical relationships as a string of 
numerical values and operators. It is useful in over- 
coming the confusion encountered with respect to mathe- 
matical operation priorities. In the example above, the 
statement : 

DO - 3.0 * A1 + 0.9 * A2 + 1 * A3 

is found. 

Confusion immediately arises as to how the sequence 
of operations should be applied. Does the statement imply 
operations grouped like: 

DO = (3.0 * Al) + (0.9 * A2) + (1 * A3), 
or is the grouping, 

DO = 3.0 * (Al + 0.9) * (A2 + 1) * A3, 
appropriate? 

Polish notation associates an operator with the two 
preceding data values. Computationally, an operation acts 
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upon the two preceeding data values in the string and the 
three are replaced in the string by the result of the 
operation. The computation begins with the first occur- 
rence of an operator in the string and terminates when 
the supply of operators is exhausted. An algorithm for 
Polish notation interpretation generally keeps the data 
values in a stack until an operator is encountered. At 
this point, the last two data items in the operator stack 
are used together with the operator to compute the result 
which is returned to the stack. Polish representation of 
the above statement might be as follows, 

DO = 3.0 A1 * 0.9 A2 * + 1 A3 *+ 
and the computation of the sequence would be 
DO = (3.0 * Al) 0.9 A2 * + 1A3 * + 

DO = (3.0 * Al) (0.9 * A2) + 1 A3 * + 

DO = ((3.0 * Al) + (0.9 * A2) ) 1 A3 * + 

DO = ((3.0 *A1) + (0.9 * A2) ) (1 * A3)+ 

DO = (((3.0 + Al) + (0.9 * A2 ) ) + (1 * A3)) 

Input to the analysis system for data transformation 
is accomplished through FORTRAN- like programming state- 
ments. These statements are restricted to arithmetic and 
logical operations on attribute, dummy and output varia-. 
bles . A syntax table (see Table IV-4) is used to define 
the statement input and when used with the syntax driver 
and a conversion subroutine with flow chart in Figure IV- 5 



TABLE IV- 4 


Data Transformation Syntax 


Transform 

Subroutine 

32 

Equation 
End of Table 

Statement 

Statement 

Equation 

Or 

Compound 


Equation 
End of Table 

Simple 

Compound 

Literal 

IF 

Subroutine 

10 


Equation 

Expression 


Subroutine 

20 


Literal 

Then 


Equation 
End of Table 

Simple 

Simple 

Equation 

Or 

Assign 


Equation 
End of Table 

End 

Assign 

Subroutine 

11 

Equation 

Variable 


Literal 

= 


Equation 
End of Table 

Expression 

End 

Literal 

End 


Subroutine 
End of Table 

12 

Expression 

Equation 

Priority 4 

Repeat 

Priority 5 


End of Table 


Priority 5 

Literal 

.EQ. 

Equation 

Priority 4 


Subroutine 

Or 

13 


Literal 

. NE« 


Equation 

Priority 4 


Subroutine 

Or 

14 


Literal 

.GT. 
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TABLE IV- 4 (cont.) 


Priority 4 


Priority 4. 1 


Priority 3 


Priority 3.1 


Priority 2 


Equation 

Subroutine 

Or 

Literal 

Equation 

Subroutine 

Or 

Literal 

Equation 

Subroutine 

Or 

Literal 
Equation 
Subroutine 
End of Table 

Equation 

Repeat 

End of Table 

Literal 

Equation 

Subroutine 

Or 

Literal 
Equation 
Subroutine 
End of Table 

Equation 

Repeat 

End of Table 

Literal 

Equation 

Subroutine 

Or 

Literal 
Equation 
Subroutine 
End of Table 

Equation 

Repeat 

End of Table 


Priority 4 

15 

.LT. 

Priority 4 

16 

.GE. 

Priority 4 

17 

. LE. 

Priority 4 

18 


Priority 3 
Priority 4.1 


Priority 3 
21 


Priority 3 
22 


Priority 2 
Priority 3.1 


* 

Priority 2 
2 3 

1 

Priority 2 
24 


Priority 1 
Priority 2.1 
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TABLE IV- 4 (cont.) 


Priority 

2.1 

Literal 

Equation 

Subroutine 

Priority 

1 

Literal 

Subroutine 

Equation 

Or 

Li teral 

Equation 

Or 

Equation 
End of Table 

Blank 


Literal 

Operand 


Equation 

Or 

Equation 

Or 

Equation 
End of Table 

Subexpression 

Literal 
Equation 
Literal 
End of Table 

Constant 


Equation 

Or 

Equation 
End of Table 

Integer 


Equation 
Repeat 
Subroutine 
End of Table 

Fraction 


Equation 

Repeat 

Literal 

Subroutine 

Equation 

Repeat 

Subroutine 

Subroutine 


* A 

Priority 1 
25 

5 

Operand 


Operand 

Operand 


" 5 ” 

Subexpress ion 

Variable 

Constant 


6 

Expression 

) 

Fraction 

Integer 


Number 

Number 

0 


Number 

Number 

0 

Number 

Number 

1 

22 



TABLE IV- 4 (cont.) 



End of Table 


Variable 

Equation 

Or 

Attribute 


Equation 

Or 

Dummy 


Equation 
End of Table 

Output 

Attribute 

Literal 

A 


Equation 

Number 


Repeat 

Number 


Subroutine 
End of Table 

2 

Dummy 

Literal 

D 


Equation 

Number 


Repeat 

Number 


Subroutine 
End of Table 

3 

Output 

Literal 

0 


Equation 

Number 


Repeat 

Number 


Subroutine 
End of Table 

4 

Number 

Subroutine 

6 


Literal 

Or 

0 


Literal 

Or 

1 


Literal 
End of Table 

9 
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Figure IV- 3. Logical Flow Diagram for SUBR.2 
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Figure IVr3. Logical Flow Diagram for SUBR.2 (cont.) 
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allows conversion of input statements to machine coded 
strings in Polish notation. An arithmetic interpreter 
subroutine is then used to effect the desired arithmetic 
operation by interpreting the Polish string prepared by 
the syntax procedures. 

Conversion Subroutine 

The conversion subroutine SUBR.2 is utilized with 
the syntax driver to convert the interpreted syntax into 
a Polish form string for subsequent computation. The 
syntax table provides the operation codes for the conver- 
sion subroutine shown in Table IV- 5 and the functions 
performed by the subroutine are listed in Table IV- 6. 

Table IV- 7 summarizes the 8-bit operation codes placed in 
the output string as a result of the syntax interpretation. 

Transformation Interpreter 

The transformation interpreter TRANS interprets the 
Polish format string provided by the IRIMAS syntax 
transformation conversion subroutine, and executes the 
indicated operations. A byte string of operation codes, 
identifiers , integer and fractional binary values and 
variable addresses constitute the input string processed 
by the interpreter. Pseudo- floating point operations are 
performed where the binary point is shifted to combine 
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Table IV- 5 

Syntax Conversion Subroutine SUBR.2 Operation Codes 


0 

Integer number 

13 

. EQ. 

1 

Fraction 

14 

.NE. 

2 

Attribute variable 

15 

. GT . 

3 

Dummy variable 

16 

.LT. 

4 

Output variable 

17 

. GE . 

5 

Constant zero 

18 

. LE. 

6 

Numerics 

19 

— 

7 


20 

STOP 

8 


21 

- 

9 


22 

+ 

10 

IF 

23 

* 

11 

ASSIGN 

24 

/ 

12 

END 

25 

** 
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Table IV- 6 

Subroutine SUBR.2 Interpretation Functions 


Operation Code 1 

The ASCII coded decimal fraction transferred 
to the LIST buffer area is converted to a 
binary fraction. The binary fraction, pre- 
ceded by a "O’* for identification, is placed 
in the output string. 

Operation Code 2, 3, 4 
These opcodes indicate than an attribute, 
dummy, or output variable has been refer- 
enced. These variables are maintained 
in arrays and the ASCII coded decimal 
number in the LIST buffer area indicates 
which member of the array is being 
referenced. The number is converted to 
a binary integer and the appropriate 
variable address is calculated. The 
variable address, prefixed with a value 
M 02 ,f to specify a variable address, is 
inserted in the output string. 

Operation Code 5 

A value of zero, prefixed with an integer 
identification ’’OO'*, is placed in the 
output string. 

Operation Code 6 

The character referenced by the input 
cursor moved to the character string at 
position LIST. 

Operation Code 7-25 

Mathematical, logical, and control func- 
tions are represented by operation codes 
7-25. These same functions are required 
by the computation software, therefore, 
the conversion subroutine simply places 
the codes in the output string. 
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Code 

0 

1 

2 

7-9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


Table IV- 7 

Polish String Operation Codes 

Function 

A 16-bit binary integer follows 
A 16-bit binary fraction follows 
A 16-bit address follows 
No operation specified 
Conditional expression (IF) 
Assign 
End 

Logical equal (.EQ.) 

Logical not-equal (.NE.) 

Greater than 
, Less than 

Greater than or equal 
Less than or equal 
No operation specified 
Stop 

Subtract 

Add 

Multiply 

Divide 

Exponentiate 
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integer and fraction values which can be handled by integer 
arithmetic operations. The specific operation codes and 
data identifiers are discussed in previous sections and 
are summarized in Tables IV- 5 (page 104) and IV- 7 (page 106)- 

Coordinate Transformation 
Transformation Definitions 

Three fundamental coordinate transformations are 
implemented by IRIMAS software. These include coordinate 
translation where a coordinate reference system is dis- 
placed along the x and y axis; coordinate rotation where 
the coordinate reference system is turned through a 
specified angle; and scaling whereby the coordinate 
reference system is expanded or contracted along the x 
and y axis. In Figures IV-4a through IV-4d are shown, by 
example, the results of these three coordinate transfor- 
mations . 

Coordinate translation is described by the following 
relationships : 

x’ = x + T x and, 

y’ = y + T’y, , 

where x* and y* are the new coordinates and x and y are 
the coordinates which are being transformed. In the 
translational mode the value of the coordinate x is 
displaced by the amount of T and the value of the 

A. 
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coordinate y is translated by the amount of T . Coordinate 

y 

transformation may be conveniently handled as matrix 
operations or vectors. For the translation case, the 
vector-matrix product becomes 


r \ 


r- * 


f * 

x' 


1 0 T 


X 



X 



y’ 

- -< 

0 1 T 

y 

► < 

y 

l 


0 0 1 


i 

< d 


i. j 


L J 


where the dummy value n l H has been introduced for 
operational convenience. Another convenience provided by 
the matrix notation arises from the understanding that a 
transformation is a single entity or operation, and that 
transformations can be combined or concatenated to yield 
a single transformation. Utilization of the matrix format 
allows ready reduction of the transformation string by 
application of the matrices on each other until a single 
matrix results. The transformation operator for trans- 
lation then is 


= < 


1 0 T 

x 

0 1 T \ 

y 

0 0 1 


The second coordinate transformation to be considered 
is that of rotation. Coordinate rotation is described by 
the relationship: 


x* = x cos 6 + y sin 0 
y f = -x sin 0 + y cos 0 
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where 0 is the angle through which the coordinate system 
is to be rotated. In matrix notation, the rotation 
transformation becomes: 



* 

cos 0 

sin 

e 

s 

0 

= i 

- sin 0 

cos 

e 

0 


0 

0 


1 


The final coordinate transformation considered is that of 
scaling. Coordinate system scaling is described by the 
relationships : 


x' = x S 

X 

y' = y s y 


where S and S are the scale factors applied to the x-axis 
x y 

and y-axis respectively. In matrix notation, the trans- 
formation is described by: 


T 


S 



0 o] 



0 1J 


Once concatenated matrices have been reduced to a 
single transform operation, the dummy variable need not be 
calculated in the final result. Thus, the transformation 
operation becomes 
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Syntax and Interpretative Subroutine 

Operator specification of coordinate transformations 
follow syntactic definitions in the same manner as data 
transformation. Data entries follow the format TR(T x , T ) 
for coordinate translation, ROT (0) for rotation and SCL 
(S x , S ) for scaling. The additional directive END 
terminates the operator entity. Examples of allowable 
operator entries are: 

TR (10,0) translate 

ROT (25) rotate 25 degrees 

TR (-10,0) translate 

SCL (10.0, 0.10) scale 

END terminate operation 

The formal syntax for coordinate transformation is defined 
in Table IV- 8. 

An interpretative subroutine, SUBR.3, is implemented 
to convert interpreted syntax into a coded string. Its 
flow diagram is in Figure IV- 9. This coded string is used 
by the TRANS subroutine (see Figure IV-6) to execute the 
coordinate transformations defined by the system operator. 
The matrix operations previously described are implemented 
in data strings within the SUBR.3 subroutine. As the 
coordinate transformations defined by the operator are 
interpreted, the appropriate transformation parameters are 
inserted into the coded strings. Subroutine operations 



TABLE IV- 8 


Syntax for Coordinate Transformation 


Transform 

Subroutine 

0 


Equation 

Or 

Operation 


Equation 
End of Table 

End 

Operation 

Equation 

Translate 


Subroutine 

Or 

7 


Equation 

Rotate 


Subroutine 

Or 

7 


Equation 

Scale 


Subroutine 
End of Table 

7 

Translate 

Literal 

TRL ( 


Equation 

Integer 


Subroutine 

1 


Literal 

9 


Equation 

Integer 


Subroutine 

2 


Literal 
End of Table 

) 

Rotate 

Literal 

ROT ( 


Equation 

Integer 


Subroutine 

3 


Literal 
End of Table 

) 

Scale 

Literal 

SCL ( 


Equation 

Integer 


Literal 



Equation 

■ Fraction 


Subroutine 

4 


Literal 



Equation 

Integer 


Literal 

« 


Equation 

Fraction 


Subroutine 

5 


Literal 
End of Table 

) 



TABLE VI- 8 (cont.) 


Integer 

Equation 

Number 


Repeat 

Number 


Subroutine 
End of Table 

9 

Fraction 

Equation 

Number 


Repeat 

Number 


Subroutine 
End of Table 

10 

Number 

Subroutine 

6 


Literal 

0 


Or 



Literal 

Or 

1 


Literal 

Or 

9 


Literal 

Or 

- 


Literal 
End of Table 

+ 

End 

Literal 

End 


Subroutine 
End of Table 

8 
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Figure IV- 5. Logical Flow Diagram for SUBR. 3 
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Figure IV- 5. Logical Flow Diagram for SUBR. 3 (cont.) 
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Figure IV-6. Logical Flow Diagram for TRANS (cont.) 
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codes are tabulated in Table IV-9. 

Data Display 

Man is accustomed to thinking in pictures and thus 
images provide him with a very personal view of his world. 
The use of thematic cartography is invaluable to the 
regional, natural resource, or land-use planner dealing 
with the spatial extent of land-use. The capability of 
image presentation is vital to the success of a resource 
information system. 

Three graphical display devices are commonly available 
for image presentation, the line printer, the digital 
plotter and a cathode ray tube (CRT) display. The line 
printer and the CRT display, for the most part, print or 
plot data on an x-y coordinate grid in a sequential manner 
rather than by graphical interpretation of coordinate 
functions. Simply, the grid or raster data format is 
most compatible with these devices. In view of this 
requirement, IRIMAS software also converts the machine 
representation of the boundary format into the grid format. 

Several functions are important to the format conver- 
sion of the boundary format data. Algorithms are available 
in IRIMAS for the tracing of boundary lines into an image 
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Table IV-9 

Syntax Conversion Subroutine SUBR. 3 Operations 


Operation Code 0 

The interpreting subroutine is initialized. 
Operation Code 1 

The value for an x-coordinate translation, 

T , is interpreted. 

Operation Code 2 

The value for a y-coordinate translation, 

T , is interpreted. 

Operation Code 3 

The value for the angle of a coordinate 
rotation in degrees, 6, is interpreted. 

Operation Code 4 

The value for an x-coordinate scale change 
is interpreted. 

Operation Code 5 

The value for a y-coordinate scale change 
is interpreted. 

Operation Code 6 

The input cursor points to a character in 
an ASCII coded decimal number which is 
moved into a buffer area for conversion. 

Operation Code 7 

The coordinate transformation specified 
has been completely translated. 

Operation Code 8 

The end of the definition phase for 
coordinate transformation is indicated. 

Operation Code 9 

The ASCII coded character string of 
numerals in the buffer storage area is 
converted to a binary integer. 

Operation Code 10 

The ASCII coded character string of 
numerals in the buffer storage area is 
converted to a binary fraction. 
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buffer area, as well as for shading the area within the 
boundary. These algorithms convert the lineal and area 
boundary format data into grid data where coordinates are 
associated with memory locations. The point- to-grid 
format conversion is direct. An additional algorithm is 
provided to clip the boundaries of the information for 
display to fit the area defined by the display extent. 

Clipping 

A clipping algorithm is implemented in IRIMAS which 
rejects those portions of the data boundaries which lie 
off of the window area defined by the display. The 
algorithm additionally computes the endpoints of the 
boundary lines at the intersection with the display boun- 
daries effectively redefining the data boundary to lie 
entirely within the display window. The implemented 
clipping algorithm, attributed to Cohan and Sutherland by 
Newman et. al. [35], not only rapidly finds the endpoints 
off the window but is also able to quickly reject any line 
which lies entirely off the display area. 

The algorithm is in two parts; the first part applies 
a rejection test. This test is implemented by extending 
the edges of the window so that the space occupied by the 
unclipped picture is divided into nine regions. Each of 
these regions has a four bit code as shown in Figure IV- 7. 
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1001 

1000 

1010 

0001 

0000 

0010 

0101 

0100 

0110 


Figure IV- 7 . Logical Codes for Clipping Algorithm 
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The endpoints of the line are assigned codes according to 
the region in which they fall. The four bits in the code 
mean the following if set: 

first bit: point is above top edge of display 

second bit: point is below bottom edge 

third bit: point is to the right 

fourth bit: point is to the left. 

If the four bit codes for both endpoints are zero, the 
line lies in the window, and if the intersection of the 
two codes is not zero, the line must lie entirely off the 
display. 

If the line cannot be eliminated by either* of these 
tests, then the second portion of the algorithm is applied. 
A line which lies partially out of the display area can be 
divided into two parts, one of which must lie entirely out- 
side the display area. A simple method of subdivision is 
to find the line's point of intersection with one of the 
window edges and to throw away that part which can be 
rejected. For example, the line AB in Figure IV- 8 could 
be subdivided at C, and applying the rejection tests the 
portion AC can be discarded. The line BC cannot be simply 
categorized and therefore must be subdivided again at D. 

The result BD is found to lie entirely inside the window 
area. In applying the rejection test repeatedly, the 
order of fhe rejections does not affect the final result. 
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Figure IV-8. 


Example of Line Rejection and Clipping 
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Line generation 

Information system data entries handled by IRIMAS are 
represented as piecewise linear boundaries specified by 
line segment endpoints. An algorithm is implemented for 
expansion of the endpoint coordinates into the set of 
coordinates which constitute the line. The algorithm is 
based on parametric representations of the x and y 
coordinates of a line. The set of parametric equations 
are expanded from endpoint to endpoint as the coordinates 
of the intermediate points on the line are calculated. 

A line with endpoints (x^ y 1 ) and (x 2 , y 2 ) can be 
represented by a set of parametric equations: 

x = + (x^ - x^)t 

y = Y X * <y 2 - ypt 

where CK t <1 . 

At the starting point of the line (x = x^., y * y^) 
t equals 0 and at the termination of the line segment 
(x = x 2> y = y 2 ) t equals 1. The line representation is 
constrained by the grid coordinate representation to have 
coordinate values at discrete integer valued points. Thus, 
the line becomes a set of integer points. This suggests 
that the parameter "t" can assume only discrete values 
corresponding to the discrete coordinates. 

In tracing the line, the algorithm generates two 
values of "t" for each additional coordinate point. The 
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value t = t 1 is computed by incrementing the x-value of 
the previous coordinate by Ax and calculating 

(x - x.) + Ax 

t = - * 

1 ( x 2 ' X P 

This computes a value of the parameter t keyed to the 
movement along the x-axis. A value for t = t 2 is computed 
by incrementing the y-value by Ay 

(y - y,) + Ay 

t = ± * 

2 (y 2 - y l> 

This value for t is keyed to movement along the y-axis. 
Since t is a monotonically increasing value from the 
initial coordinate (x^ , y^) the coordinate point of the 
line closest to the preceding point has the smaller t 
value. Thus, the next point in the coordinate sequence 
tracing the line is generated by the smaller of [t^ , t 2 J. 

Since the coordinate sets (x^, y^) and (x 2 , y 2 ) are 
arbitrary, movement from x-^ to x 2 , and y^ to y 2 can require 
either positive or negative increments depending upon the 
relative relationships of the coordinates. Thus, the 
value of Ax and Ay used to increment the line coordinates 
can either be positive or negative depending upon the 
particular direction of the line. Appropriately chosen, 
the parameter t remains positive valued between 0 and 1. 
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Area Generation 

The display of area information requires that the 
coordinate values of points within data boundaries be 
identified and associated with the information values of 
the area. In IRIMAS, an algorithm is provided to determine 
those points internal to a boundary and to compute the 
intersection of these points with the display window. 

This '’shading" algorithm computes the intersection of 
the line y - y^, with the area boundary segments for values 
of y^ within the upper and lower y-limits of the line 
endpoints. These intersections are compared with each line 
segment's endpoints to determine if the intersection is on 
the boundary of the area. These x-coordinate values are 
sorted in ascending order and are compared with the x-axis 
limits of the display area. Those points without the 
display area are replaced by the display limits. The 
resultant set of x-axis values, taken in pairs, provide 
upper and lower limits of the points along the line y = y^ 
which are both in the bounded area as well as in the dis- 
play window. In Figure IV-9 is shown a graphic example 
of this shading performed by the algorithm. This algorithm 
provides a method for identifying the points encompassed 
by the data entry area which lie with the display area in 
order that area information can be transformed into the 
grid format. 
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Figure IV- 9. Example of Algorithmic Shading 
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Expansion Subroutine 

An IRIMAS expansion subroutine EXPAND is implemented 
to convert the machine format data into the grid format for 
display and image- function operations. This subroutine 
implements the described algorithms for clipping, line 
generation and area generation and logically intersects 
the generated coordinate values with specified coordinate 
bounds. If a calculated data point lies within the 
specified coordinate frame, the subroutine execution is 
interrupted. At this point, the data values may be 
inserted into a coordinate indexed array thus completing 
the expansion. In Figure IV- 10 is shown a flow graph for 
EXPAND. 

Optionally, the subroutine interruption may be 
utilized as an indication of congruence between the machine 
format data set and an additional data set used to define 
the data frame input to EXPAND. For example, rectangular 
classification training sets defined as data frame boun- 
daries cause an interrupt to occur when the input data 
corresponds to the training set areas. Under this con- 
dition, data statistics could be computed for utilization 
in classification algorithms. Alternatively, coordinate 
indexed locations in an array could be tested against a 
preset logical condition to control further processing. 

The subroutine EXPAND provides for convenient conversion 
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Figure IV- 10. Logical Flow Diagram for EXPAND 
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Figure IV-10. Logical Flow Diagram for EXPAND (cont.) 





Figure IV- 10. Logical Flow Diagram for EXPAND (cont.) 
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of the machine representation of boundary data to coordi- 
nate indexed arrays for further image processing or 
display . 
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CHAPTER V 
IMAGE PROCESSING 

The current availability of timely synoptic imagery 
from large portions of the productive land and ocean areas 
of the world offers an unprecedented source of information 
to the resource manager. Properly interpreted, this 
information can provide new insight into resource avail- 
ability and productivity and especially the environmental 
impact of human interaction with natural ecosystems. In 
the preceding chapter, an information storage and retrieval 
format is defined to provide for effective and flexible 
manipulation of data related on a spatio-temporal basis. 
Algorithmic procedures are described for the manipulation 
and analysis of these aggregated data. A revised method- 
ology is required for the development of these information 
bases from the newly available earth imagery. 

Large quantities of significant resource information 
are embedded in earth imagery, which is available on a con- 
tinuing basis, particularly through the NASA Earth Resources 
Technology satellite systems. Consequently, major image 
analysis and interpretation aids are implemented in IRIMAS 
to enable an interpreter to interactively enhance, classify, 
transform and display the image product. Additional soft- 
ware is provided to allow the analyst to save his interpre- 
tation in the format of the IRIMAS data files. These major 
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image analysis components provide the capability for the 
entry of image interpreted resource information into the 
complete IRIMAS system. 

Image Analysis Functions 

The implemented image analysis functions are ones 
which have been proven by repeated usage for resource 
survey and reconnaissance applications. Image enhancement 
by smoothing and edge enhancement techniques are estab- 
lished as useful procedures in the detection and delinea- 
tion of specific image features [36]. The smoothing 
functions reduce the effects of random image noise on the 
interpre tability of the data, while the edge enhancement 
functions provide a method whereby boundaries inherent in 
the image are enhanced and their detectability improved. 
These image enhancement techniques are particularly 
restricted to single spectral bands in imagery and do not 
offer the flexibility for the handling of multiband 
information (extensions to these multidimensional data 
sets, however, seem reasonable). 

Multispectral imagery requires a more complicated set 
of procedures, specifically clustering or classification. 
These analysis techniques essentially compute similarity 
functions between data points and group data points which 
are similar into categories which are dissimilar from 
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other categories. The similarity functions may be either 
deterministically or stochastically derived. 

An additional enhancement technique, density slicing, 
is beneficial in those cases where the magnitude resolution 
of the data is restricted either by system or interpre- 
t ability constraints. This technique requires the ordering 
of the data into groups which contain equal numbers of 
elements. Assignment of values or symbols to each group 
completes the ’’density slicing” transformation and improves 
the ability to distinguish homogeneous categories by 
reducing the confusion generated by interpreting image 
elements with a wide range of values. 

Smoothing 

The purpose of the image smoothing function is to 
suppress or average the random noise introduced into the 
data by the measurement system and to reject and replace 
those image elements which appear to be improper with 
respect to the surrounding picture elements. The smoothing 
operation which is implemented is called noise-point elim- 
ination [4] . This is a "neighborhood" operation where 
several points in the vicinity of the element being tested 
are used in the smoothing operation. The noise-point 
elimination technique replaces a point by the average of 
its neighbors if the point value exceeds that average by 
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some preset threshold. This technique has the least 
effect on class boundaries while still providing some 
correction for image data noise. Additional approaches to 
the smoothing operation include replacing a specified 
point by the average of its neighbors, however, this 
approach tends to soften the edges in an image and thus 
tends to be detrimental to resource interpretation. 

Edge and Line Enhancement 


The enhancement of edges, lines and boundaries in an 
image is significant for improved interpretation. Broad 
categories can be brought into sharper distinction through 
application of line and edge enhancement algorithms. An 
item of caution in edge enhancement concerns the enhance- 


ment of lines and edges with regard to a particular 
direction. The results of some enhancement procedures are 
affected by the particular orientation of the edge fea- 
tures. A suitable enhancement algorithm which reduces 
the effect of this directional dependence is the gradient 
of the image function at the point being enhanced. In 
this operation, the gradient 


f(x,y) - 


(x,y,) 2 + 3f fx,y) 


1/2 


3x 9y 

is computed to replace the data point. Computation of 
this function for digital pictures can be approximated by 



where is the i-th picture element from the j-th row. 

The gradient operation, and for that matter differen- 
tiation or differencing as a whole, directly relates to 
the concept of spatial filtering prevalent in image 
enhancement literature today [36], Two dimensional spatial 
frequencing domain transformations, such as the Fourier 
transform, are used to convert the image function from the 
spatial domain to the spatial frequency domain. Selection 
of the high frequency components in the spatial frequency 
domain and inverse transformation emphasizes or enhances 
those areas in the image where the image function is 
rapidly changing as a function of the coordinate variables. 
Thus, the enhancement of the image occurs at the boundaries 
or edges between categories of varying image densities. In 
an analogous manner, smoothing effectively selects and 
enhances the low spatial frequency components of the image, 
essentially providing low-pass filtering. 

Density Slicing 

The two-variable (dimension) image function which 
assigns a reflectance value to each coordinate point in 
the image plane, typically is quantized into numerous 
levels. The ERTS-1 MSS system for example quantized its 
image measurements into 64 distinct levels. The vari- 
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ability in scene radiance, signal to noise ratios, and 
quantization levels all affect the distribution of 
picture elements in the range of this discrete function. 
Manual image interpretation cannot effectively utilize 
this wide range of information values. Consequently, a 
transformation of the fine density resolution into a more 
suitable (coarser) range for manual image interpretation 
is appropriate. The most obvious procedure for effecting 
this transformation is to linearly modify the set quan- 
tization levels such that a smaller number of density 
values (e.g. 8) are produced for interpretation purposes. 

This linear transformation is unacceptable due to the 
wide variablility in the scene illumination and also in the 
reflectance of specific categories of interest. Typically, 
the image function for a specific scene may only cover 
20-303 of the allowable quantization range. A linear 
transformation then assigns that percentage of the 
available quantization levels to the whole image, 
resulting in a significant loss of image information. A 
non-linear transformation seems appropriate. 

The image density transformation which provides the 
greatest preservation of information, non-linearly sets 
new quantization levels according to the distribution 
of the image points in the range of the discrete image 
functions. These assignments are made so that equal 
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numbers of image points are assigned to each quantiza 
tion level. Constructing the quantization levels in this 
manner allows each image density value to have an equal 
weighting with respect to every other image point. 

Clustering 

Clustering is normally thought of in terms of 
classification, since in general this category of tech- 
niques groups image points into classes where all the class 
members have similar image density values and non-members 
of a class are dissimilar. Numerous clustering algorithms 
have been described [6] and utilized as classification 
tools. The major drawback to the use of clustering as a 
classification mechanism is the large amount of processing 
required for multiple passes through the data in search 
of some suitable collection of categories with adequate 
classification accuracy. Additionally, the classes 
selected by a clustering algorithm generally do not assure 
that the classes determined solely from the spectral 
distribution of measured reflectance are those of concern 
to a resource manager. Consequently, the clustering 
techniques have been to a large extent shadowed by other 
classification techniques, for example, maximum likelihood. 
Clustering does offer specific advantages where the 
accumulation of similar points and the assignment of a 
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single grey scale to these groupings may be useful as an 
enhancement aid. IRIMAS uses a single pass clustering 
algorithm to group similar image points which lie together 
in a general area. 

The clustering technique implemented [37] provides a 
reasonable trade-off between classification and speed of 
computation. The algorithm also uses a norm for the 
computation of point similarity which is more reasonable 
in terms of the actual classification assumptions than the 
maximum likelihood or a similarly derived ’’least squares” 
categorization technique. The distance function used for 
the clustering algorithm is the & 00 metric d, 

d - max | x^ - | 

where the distance in the measurement space is the maximum 
of the absolute values of the differences between the 
spectral components. It was found in [37] that this 
distance measure provided an increased sensitivity to 
spectral change, suppressed the smearing of categories in 
the measurement space, and its success could be more 
readily predicted from individual data band characteristics 
than other measures averaging all the bands. 

Spatial clustering is determined by computing the 
spectral cluster center and by adding points to the cluster 
where the distance from the point to the cluster center 
was less than a specified threshold distance away. The 
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clustering algorithm with its associated measure is imple- 
mented in XRIMAS to provide local enhancement of class 
associations . 

Many resource categories of interest have elements 
which are aggregated spatially. The implemented clustering 
algorithm builds a table of clusters as it processes 
through the data. A new data point is compared to the 
cluster to which the previous point belonged, if the new 
data point is determined to not belong to that cluster 
the point is compared to the cluster of the point that 
lies near it on the preceding scan line. For typical 
cluster- class if ication, if this test fails, the point 
would be compared to all existing clusters before a new 
cluster was initiated. However, for the enhancement mode 
in which this procedure is utilized, the maintenance of 
a large number of clusters is generally not desirable. 
Therefore, if the point does not belong to the clusters of 
either of the two preceding points, a new cluster is 
initiated. At the end of each scan line, the clusters 
not used for that scan line may be reinitialized for use 
since the point aggregations which defined them are 
disjoint from the subsequent data points. 

This cluster enhancement algorithm has the net effect 
of smoothing data within spatial aggregation while 
enhancing their boundaries with other groupings. Such a 
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procedure is helpful in identifying major data groupings 
particularly to aid in the selection of training sets for 
a supervised classification mode. 

Classification 

A supervised classification procedure is utilized by 
IRIMAS to provide class selection and identification from 
mult ispectral image function inputs. In a supervised 
classification mode, the classes to be identified are pre- 
specified to the algorithm which compares the new data to 
these specified classes and determines which class best 
describes the point. The classification procedure which 
has produced the best classification results operationally 
is the maximum likelihood classification procedure applied 
at the Laboratory for Applications of Remote Sensing, at 
Purdue University [28]. This algorithm has been generally 
accepted and versions of the LARSYS series of software 
packages are in wide usage. Consequently, the maximum 
likelihood classification procedure is implemented in 
IRIMAS. 

The maximum likelihood classification algorithm 
assumes that the data fall into categories which may be 
described by multivariate Gaussian distributions with 
differing means and covariances. The implementation of 
the classifier algorithm computes the value of the con- 
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ditional probability density function for each category 
to be considered, and infers the membership of the data 
point to that category for which the conditional density 
function is the largest. 

Class means and covariances are estimated from the 
measurement of signatures which are known or assumed to 
belong to that class. These estimated statistical 
parameters are then used in the classification algorithms. 
A class mean (for class i) is defined to be 



n 

and the estimated class covariance 


is 


. E . = -A- E(X. “ X.) (X. - X.) T where 
1 _ 1 i J K l i J 

n- 1 

X^ is the n~ dimensional measurement signature for the i-th 
class. The conditional probability then becomes 

f(X|<o.) = §exp [-1/2 (X-X.) T e'V-X.)] 

1 (2ir lEj) 2 1 1 

where X is the measurement signature, is the i-th class 

and X., and Z. are the estimated i-th class means and 

i* l 

T 

covariances. The matrix operations ( ) , | | and 

( ) ^ correspond to transpose, the determinant, and the 

inverse of a given matrix. The classification decision 
for two classes then is 


to 


to Q if f Cxjto Q )> : f (xlto^) 
to i if f (x|aO<f (x|u> o ) 
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for the data point X. An additional test is inserted into 
the classification algorithm. The possibility exists that 
a data point does not belong to any of the data classes 
specified. In this case, the data point is expected to 
have a low probability of belonging to any of the classes. 
This point is then assigned to the null or zero class 
where the largest value of the density function is 
determined to fall below a preset threshold. Using this 
rejection class minimizes the chance of identifying a 
measurement as belonging to a class when in fact the class 
to which it belongs was not specified. 

Image Processing Software 

A software component entitled PROCES is implemented 
in IRIMAS to provide the image processing functions or to 
reference specific subroutines supporting application of 
these analysis techniques. This software driver has two 
specific functions, the conversion of data byte strings 
into integer arrays for processing; and to effect the 
analysis function through either calls to external sub- 
routines, or internal implementation of the analysis 
techniques . 

The operation of the software driver PROCES is 
controlled by basically two parameters. The first param- 
eter PASS indicates whether the driver should be in an 
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initialization, processing or termination mode. Entry 
into PROCES with PASS equal to 1 causes the driver to 
initialize and zero the appropriate arrays and variables. 
Entry with PASS equal to 0 indicates normal processing of 
the input data values, and entry with PASS equal -1 in- 
dicates that processing is concluded and termination of 
the analysis sequence is to be undertaken. The second 
control parameter identifies the specific analysis tech- 
nique to be applied. In Table V-l, is shown the operation 
codes and the specific enhancement or classification 
technique to be applied. 

In addition to the processing techniques described, 
the IRIMAS software system allows operator definition of a 
special transform to be applied to the data. This 
implementation utilizes the syntax drivers and interpre- 
tative subroutines previously described in Chapter IV for 
specification of data transformations. The capability 
provides a means for the application of specific discipline 
related transformation without major system modification. 

An operational flow chart for PROCES is shown in Figure 


V-l. 
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TABLE V-l 

PROCES Operation Codes for Enhancement and Classification 


Code 

0 

1 

2 

3 

4 

5 

6 
7 


Function 

No Operation 

Classifier Training 

Maximum Likelihood Classification 

Algorithmic Clustering 

Edge Enhancement 

Smoothing 

Histogram Preparation 
Operator-defined Transformation 





^ RETURN J 

Figure V-l. Logical Flow Diagram for PROCES 
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Figure V-l. Logical Flow Diagram for PROCES (cont.) 






Figure V-l. Logical Flow Diagram for PROCES (cont.) 









Figure V-l. Logical Flow Diagram for PROCES (cont.) 
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Figure V-l. 


Logical Flow Diagram for PROCES (cont.) 






Figure V-l. Logical Flow Diagram for PROCES (cont.) 
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CHAPTER VI 
SYSTEM OPERATION 

The preceding chapters discuss the general require- 
ments for a resource information management and analysis 
system to develop an effective interface between a resource 
planning function and the information available from 
remotely sensed data. An information format and techniques 
for information retrieved from data stored in that format 
are developed. Image analysis functions which enable an 
interpreter to analyze satellite imagery interpreting 
resource information, and to save it in data files for 
future reference and analysis are also discussed. Addi- 
tionally, the software subroutine and techniques developed 
for these analyses and information data handling are 
described. This chapter concludes the detailed description 
of data handling and analysis techniques implemented, by 
explaining the integration of these individual subprograms 
into the Interactive Resource Information Management and 
Analysis System, IRIMAS, and by further defining the actual 
interface with the analyst/system operator. 

System Control 

The software system IRIMAS has three principal 
function areas: 1) input; 2) processing, and; 3) output, 

which are operator specified. The input function takes 
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data from an input data file and converts them into a 
usable internal format for IRIMAS. The processing function 
implements data and coordinate transformation, as well as 
image enhancement and classification algorithms. The out- 
put function provides image display, printer tabulation, 
and information format output to IRIMAS data files. 

Operator interface has been accomplished in two pri- 
mary ways. Operator directives are provided to the system 
for definition of the system tasks and are interpreted 
according to a syntax table in a manner similar to the 
procedures previously described for syntax analysis. The 
interpreted syntax provides operation codes which are saved 
for further interpretation by the system to define specific 
system operations related to data input and analysis. 
Specification of the output data product, or data display 
is also provided. . Operation directives for control and 
termination of the system tasks are also implemented. 

During the course of system analysis operations, 
operator inputs are also required to define analysis 
parameters. These are obtained by computer initiated 
queries to the system operator. The response to these 
questions provides the required information for the system 
analyses to continue. This approach to communication 
allows the computer to generate sufficient descriptive 
narrative accompanying each request to familiarize the 
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operator with the data input requested. 

Operator directives for system task specifications 
define input, processing, output and control functions. 

These directives are saved in a stack array and are 
sequentially executed by the system. Additionally, the 
processing functions of IRIMAS are controlled by a loop 
parameter, PASS. This parameter specifies whether the 
system is in an initialization mode prior to the input 
of data (PASS=1) , an operation or processing mode (PASS=0) 
where the data is available for processing, or a termination 
mode (PASS=- 1) for the completion of processing at the 
end of the available data. The pass parameter is primarily 
controlled by the input functions which set the parameter 
according to the availability of data to the processing 
subsystem. 

Analysis processing in IRIMAS is accomplished by each 
system element on single data entries which are provided 
by each call to an input driver. Consequently, a loop 
directive, "GO TO is installed to set up the processing 

loop back to the driving input function. The loop is 
broken when the input function detects an end-of-data 
status and sets the PASS parameter to the termination mode. 
At this time, operations subsequent to the input function 
are executed in the termination mode, however, the control 
branch does not occur in the termination mode. Operator 
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control directives PAUSE and STOP halt the execution of 
the specified functions. The major difference between 
these functions is that the system must be restarted after 
a STOP, and for a PAUSE the system resets to input addi- 
tional operator function directives ^ The allowable opera- 
tor directives are tabulated in Table VI- 1, and are 
explained in following paragraphs. 

The system IRIMAS is developed on a minicomputer 
system configured with printer and color CRT output 
displays. Restricted by the application of the display 
techniques, the system generates an output format with 
both an image and alphanumeric representation. It is 
anticipated that in general, the color CRT system is not 
readily available to all minicomputer facilities and 
therefore the output products have been designed for 
display on line printers. Several output display formats 
are provided by IRIMAS. These are primarily graphics or 
image output to either the printer or to the CRT, and 
tabulations of analyses parameter providing on the printer. 

IRIMAS Input Drivers 

The IRIMAS system accepts two basic varieties of input 
data. The first input variety is the standard IRIMAS 
information data format described in Chapter IV. In this 
format, data are saved as strings of characters specifying 
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Directive 

INPUT 

ERTS 

TRAIN 

CLASSIFY 

CLUSTER 

EDGE 

SMOOTH 

SLICE 

TRANSFORM 

XDATA 

XCOORD 

EXPAND 

DISPLAY 

MAP 

DEFINE 

HISTOGRAM 

STATISTICS 

DENSITY 

OUTPUT 

STOP 


TABLE VI- 1 

IRIMAS Operator Directives 
Function 

Input from IRIMAS files 

Input from ERTS tapes 

Classification Training 

Maximum Likelihood Classification 

Algorithmic Clustering 

Edge Enhancement 

Smoothing 

Histogram Preparation 

Interactive Image Data Transformation 

Interactive Information Data Transforma- 
tion 

Interactive Coordinate Transformation 

Algorithmic Expansion of Information 

Output to CRT 

Output to Printer 

Set up IRIMAS Files 

Print Histogram 

Print Classification Parameters 
Print Density Slice Levels 
Print Output Variables 
Halt IRIMAS Processing 
Pause in IRIMAS Processing 


PAUSE 
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GO TO # 


TABLE VI- 1 (cont.) 

Transfer to Statement Indicated 
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resource boundaries by coordinate locations, attribute 
values, and data set descriptions. It is the format which 
allows the interpretation of images to be uniquely saved 
and referenced without the intermediate interpretation 
steps. This format is also convenient for the inclusion 
of additional resource information obtained from ground 
survey or other techniques, provided the information can 
be referenced to a common coordinate basis. The software 
driver INPUT implements this function. 

The second input source supported by IRIMAS is from 
the Earth Resources Technology Satellite. This data input 
is from digital magnetic tape and is in a fixed raster 
scan format. It contains the four spectral bands of 
reflectance data measured by the ERTS satellite multispec- 
tral scanner system. The IRIMAS component "ERTS” provides 
the software support for the reading of ERTS digital data, 
and converting the ERTS format into an internal format for 
processing . 

INPUT 

The software driver INPUT implemented in IRIMAS 
facilitates resource data input from IRIMAS data files. 

The information format is defined by a syntax table, and 
syntax interpretative techniques are utilized to convert 
information from the data files into a machine internal 
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format. Fundamental to the system is the internal buffer 
area where the information from the input data file is 
kept. This buffer area is named INFORM and includes three 
specific arrays. The first array, named ATTR, holds the 
values of the attributes which refer to the aggregate of 
information elements being processed from the input file; 
the second array, COORD, holds the x-y coordinate boundary 
of the category; and the array NAME holds the name of the 
data entry. This internal format is fixed to allow effec- 
tive processing of the referenced information by processing 
subroutines . 

ERTS 


Input to the IRIMAS processing system from the Earth 
Resources Technology Satellite digital image tapes is 
accomplished by the software driver ERTS. This subprogram 
allows operator selection of the tape line and data cell 
coordinates of the ERTS tape and also the selection of 
either one of the four ERTS bands from the tape or all of 
the bands. After the selection of the ERTS input param- 
eters, the subprogram reads in the ERTS data, a line at a 
time, converts the data according to the bands selected 
and stores the data into a raster buffer area named BUFFER. 
This buffer array holds two lines of ERTS data simultane- 
ously to allow neighborhood enhancement processing. 



164 


On the initial entry to the input driver, ERTS tape 
line, data cell and band parameter specifications are 
requested from the operator. The subprogram keeps a count 
of the number of tape lines read from the tape and either 
skips or backspaces to the line requested before processing 
begins. Successive calls to ERTS cause a fetch of an 
additional data raster line and provides for processing 
into the buffer area. After 128 lines of ERTS data are 
processed, the driver interrupts the processing loop with 
an end-of-data indication which sets the PASS parameter 
for the termination mode. 

Processing 

The second major set of functions implemented in 
IRIMAS are concerned with the processing of the input data. 
The processing drivers deal with one of three data storage 
arrays. The first two arrays, INFORM and BUFFER, were 
previously described (Chapters IV and V) as internal format 
storage for data input from IRIMAS data files and from 
Earth Resources Technology Satellite digital tapes, 
respectively. The third array area, FRAME, is used as an 
intermediate storage buffer for grid format data and 
corresponds to an output display window. Information 
inputs from IRIMAS data files may be expanded into a grid 
format and saved in this storage array. Different infor- 
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mation entries may be overlapped in the array and analysis 
performed. Additionally, the expansion into the frame 
array provides a convenient format from which data may be 
displayed. 

Processing of IRIMAS input data is interactively 
specified by an operator. The input commands to the system 
together with their meanings and referenced processing 
drivers are described in subsequent paragraphs. 

TRAIN 

The TRAIN operation allows the computation of mean 
value and covariance matrix estimates which are used in the 
maximum likelihood algorithm. Initialization procedures 
for this operation request information from the system 
operator to set the number of classes to be specified, 
and the names and rectangular boundary coordinates for 
each class. The coordinate of each data point during the 
processing mode are compared to the coordinates specified 
for the class boundaries. For those points lying within 
the boundaries, the values for each of band specified by 
ERTS input processing are accumulated into buffers in 
preparation for the application of the statistics estima- 
tion formula. Entry into this software processing routine 
during the termination pass causes the statistics estimates 
to be completed and the class statistics to be saved. 
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CLASSIFY 

This classification directive specifies maximum like- 
lihood classification of the input data. During the in- 
itialization mode, the covariance statistics computed by 
the TRAIN processor are inverted and their determinant 
calculated for classification processing. During the 
process mode, each input data point is utilized by the 
classification processor to determine its most likely 
class and to compare the probability of its being produced 
by that class within a preset threshold. If the probabil- 
ity of the point having occurred, (conditioned on the most 
likely class) , exceeds the threshold; the first point 
attribute, or band value, is set equal to the class number, 
otherwise, it is set equal to zero. The number of the 
bands utilized in the. classification mode is determined by 
the bands specified in the ERTS input driver. Termination 
processing has no effect on the classification procedure. 

CLUSTER 

The system directive CLUSTER provides for the 
clustering of adjacent picture elements whose computed 
similarity is below a preset threshold. Initialization 
causes each set of cluster parameters to be zeroed. Note 
that the cluster parameters are saved in the same buffer 
area as the training set boundaries so that the two 
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processing techniques may not be applied simultaneously. 
The processing mode of the cluster algorithm computes the 
clusters to which incoming data vectors belong and output 
the cluster number in the position of the first input data 
band. The termination mode does not affect the cluster 
algorithm and the number input data bands processed is set 
by the ERTS routine. 

EDGE 

The system directive EDGE specifies edge enhancement 
of the input data according to the gradient technique 
described in Chapter V. A neighborhood of three is 
utilized, consisting of the data point and adjacent 
neighbors which occur later in the input data string. 
Single band data are assumed and the edge enhancement 
operation is not affected by initializing or termination 
modes . 

SMOOTH 

The system directive for smoothing provides for 
noise-point elimination by comparing a data point with the 
average of its four-point neighborhood. If the data point 
is significantly different, exceeding a preset threshold, 
the data point is replaced by the neighborhood average. 

A single band is assumed, and initialization and termina- 
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tion modes are ignored. 
SLICE 


The computation of a 128 value- level histogram is 
specified by the SLICE directive. The break points for 
density slicing the data values into equal number bins are 
also computed. Initialization causes the histogram data 
arrays to be zeroed. During the processing mode, the 
array elements are incremented according to the value of 
the data point processed. Only one data band is assumed. 
During the termination mode, the density slice breakpoint 
levels are calculated. 

TRANSFORM 

The system directive TRANSFORM sets up interactive 
programming of input data transformations. Upon initial- 
ization, the processor requests a definition of the trans- 
forms to be utilized. The input statements are checked 
against the data transformation syntax described by SNTBL2 
and interpreted by SUBR. 2 previously described. The 
transformation definition terminates with the input of the 
"END" statement. The interpretation of the defined trans- 
form is saved and during the processing mode is utilized 
by subroutine TRANS to effect the data transformation. 

The number of bands available for processing depends upon 



the specification provided by ERTS. The bands are 
referenced as attribute variables (Al, A2 , etc.)* The 
termination mode is ignored by the transformation proces- 
sor . 

XDATA 

This operation specified by XDATA provides for the 
transformation of data input from IRIMAS data files into 
the INFORM buffer arrays. During the initialization 
procedure, the driver calls for operator specification of 
the data transformations according to the data transfor- 
mation syntax defined by SNTBL2 and interpreted by the 
subroutine SUBR.2. During the processing mode, the input 
data are transformed in accordance with the defined trans- 
formation string interpreted by the subroutine TRANS. 

The termination mode does not affect the XDATA driver. 

XCOORD 

The processing directive, XCOORD, provides for the 
coordinate transformation of data input from the IRIMAS 
data files into the storage array INFORM. Upon initial- 
ization, the driver requests the coordinate transformation 
definition from the operator and sets up the transformation 
string to be interpreted by TRANS during the processing 
phase. During processing, the coordinates are unstacked 
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from the INFORM coordinate stack and are transformed by 
the interpreter TRANS. Upon completion, the transformed 
coordinates are returned to the INFORM coordinate stack. 

The termination mode does affect the XCOORD processor. 

EXPAND 

This processing operation specifies the key routine 
in the processing of IRIMAS data files. The EXPAND 
software implements the clipping, line generation, and 
area generation algorithms described in Chapter IV for 
the conversion of IRIMAS information files into grid 
format data. These grid format data are saved in the 
storage array FRAME, and interfaced directly with the 
display peripherals. The processor EXPAND ignores the 
initialization and termination modes. 

System Output 

The final set of major functions implemented in IRIMAS 
provide for the output and display of the processed data. 
Output processing components together with the numonics 
utilized for operator specifications are described, 

DISPLAY 

The output directive DISPLAY provides for the CRT 
(cathode ray tube) display image data lines maintained 
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in the BUFFER array. These image lines are displayed on 
the system dynamic color display [38], currently in a 128 
x 128 picture element array. This output processor is 
typically specified within an input-processing loop since 
only single lines of data are output. Neither initiali- 
zation nor termination processing modes affect this 
processor. 

MAP 

The output directive MAP provides for the output of 
image data to a printer device. Density slicing of the 
data and assignment to alphanumeric symbols according to 
parameters provided by the SLICE processing provide a 
hard copy of the data image. This processor nominally 
provides a 64- column image display however, larger line 
lengths are printed in multiple parts. The processor is 
unaffected by the initialization or termination modes. 

DEFINE 

The output element DEFINE is a key module implemented 
in IRIMAS and provides for the interactive specification 
and output of IRIMAS data files. The primary elements for 
the construction of an IRIMAS data file are requested from 
the operators. These data are imbedded into the sentences 
prescribed by the data file syntax and outputs these 
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sentences as records to an output file. This module does 
not recognize any of the processing modes. 

HISTOGRAM, STATISTICS, DENSITY, OUTPUT 

Implemented output modules provide printed hard copy 
output from the analysis operations described. The 
directive HISTOGRAM, and DENSITY specify the printing of 
a graphic display of the data histogram compiled by the 
processing module SLICE, and for the output of the break 
points computed by that processing function. The directive 
STATISTICS specifies output of the data means and co- 
variances computed for the data classes defined in TRAIN, 
while the directive OUTPUT calls for the printing of the 
data contained in the output data arrays referenced by 
operator defined data transforms. The referenced output 
routines do not recognize the different processing modes. 

Control 

Three control directives are provided for operator 
control of processing. These directives include STOP, 

PAUSE and GO TO #. The STOP directive causes the program 
to terminate when this directive is encountered. This 
requires restarting of the IRIMAS system, which initializes 
certain processing variables, and zeros data arrays. The 
restart of IRIMAS also rewinds data input and output files. 
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For a temporary interruption in data processing, the 
operator directive PAUSE should be specified. This 
directive causes termination of the processing sequence 
when it is encountered, however, the interruption provide 
for a start operator definition upon system processing 
directives, as opposed to the complete restart required 
by STOP. 

Since the input drivers and the processing modules 
handle a single data record at a time, the operator direc- 
tive GO TO # has been established to set up the loop. This 
directive is used to loop back to the statement number 

generally the input directive for multiple record 
processing. The GO TO # driver resets the processing mode 
from "initialization" to "processing" and branches for 
these two modes, however it is processing mode 0, a no 
operation function in the termination mode, which is set 
on completion of the input of data required. 

Chapter VII describes by example, the operation and 
operator interface of the software system IRIMAS. Operator 
information statements and diagnostic messages are also 
discussed. 
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CHAPTER VII 
SYSTEM DEMONSTRATION 

The general specifications and specific implementa- 
tion of the Interactive Resource Information Management 
and Analysis System, IRIMAS, were discussed in detail in 
previous chapters. The system is implemented in two parts, 
one dealing with data stored in the IRIMAS format, and the 
second dealing with the image analysis functions. This 
dual implementation is required by the memory size con- 
straints of the current' configuration of the Texas Instru- 
ments 980A minicomputer system on which the software was 
developed. This chapter shows by demonstrated example 
the features of the system, the interactive operator 
interface, and the diagnostic and error messages implemen- 
ted. 

The operator interface previously described is 
implemented by computer- ini tiated dialogue and operator 
response. The system requires operator definition of 
analysis functions and input parameters during phases of 
specific function execution. Diagnostic phrases and 
system error messages are output as required to provide 
operator assurance of orderly data processing and to point 
out procedural errors in processing. 

The system is demonstrated in a manner consistent 
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with the two-part implementation. One demonstration deals 
with the information format processing and a second 
demonstration with the image processing function. A third 
demonstration on the use of the IRIMAS system in change 
monitoring is also included. 

Information Processing 

The portion of IRIMAS implemented to handle and 
analyze data in the IRIMAS data format is called INFORM. 
This program segment contains the syntax drivers for 
interpretation of data file information, and data and 
coordinate transformations. Also included is the software 
described in Chapter IV required to expand the boundary 
format representation into grid format representation. 

The output functions for CRT display and line printer 
mapping are also included. Data input is restricted to 
IRIMAS format data files and the software for creating 
these files is also implemented. This segment of the 
demonstration is intended to show the generation of an 
IRIMAS data file describing the well-known block TAM 
insignia of Texas A§M University. This insignia is used 
to represent resource information which might be specified 
in terms of the IRIMAS data files. It is recognized that 
this information relates to either points, lines, or areas. 
Subsequent productions of the IRIMAS system demonstrate 
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the reconstruction and display of the "resource informa- 
tion", as well as coordinate transformations which may be 
applied so that information from different sources may be 
aligned. Resource information such as the major highway 
systems in Brazos County may be similarly entered into the 
system and analyzed. Also shown in this demonstration are 
the operator directive sequences required for specification 
of these analyses. 

Operator Directives 

IRIMAS operator directives are required to specify 
the system functions. In Figure VI I- 1 is shown an example 
of the computer dialogue for several of the control direc- 
tives. A syntax interpreter is utilized by the software 
system to implement system function directives. 

The first entry demonstrated is an incorrect directive 
to the system. A syntax error is detected and the operator 
is instructed to re-enter the numbered line. The correct- 
ed line is entered and the software steps to the second 
line for directive entry. 

In the second line is entered the symbol "@" which 
instructs the system that the operator desires to re-start 
the function definition phase. At this point the system 
resets, and re-initiates the definition phase. 

The PAUSE directive and the STOP signify termination 
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Figure VII-1. Example of Operator Dialogue 
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of the function definition phase and indicate the function 
sequence termination point. The difference in the two 
directives is that STOP causes termination of the system 
execution. The directive PAUSE signifies that the func- 
tion definition is completed, however, the data contained 
in system buffer areas should be retained as further 
processing may be required. The system, upon completion 
of a PAUSE directive, indicates that the defined functions 
have been completed and that additional function definition 
is required. 

Completion of a STOP directive causes an operator 
diagnostic to be printed to indicate normal system 
termination and halts the IRIMAS execution. These control 
functions are also demonstrated in Figure VII- 1. 

Data File Definition 

The software subsystem DEFINE is developed to allow 
interactive operator definition of IRIMAS data files. 

This driver requests data information from the operator, 
formats the information into that specified by IRIMAS and 
data files. This driver requests data information from the 
operator, formats the information into that specified by 
IRIMAS and writes the data to an output file. 

In Figure VII-2 is demonstrated the procedure for 
definition of IRIMAS data files. In this case, the opera- 
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VII- 2. Procedure for Definition 
of IRIMAS Data Files 



tor (Step 1) specifies the system function DEFINE, and 
conclusion of the system analysis function definition by 
the PAUSE statement. 

Execution of the DEFINE function prints the explica- 
tive message which states that information data set defi- 
nition is being accomplished (Step 2). Additionally, it is 
explained that individual data entries are terminated by 
the entry of a blank line, and the entire data set may be 
terminated by an end-of-file symbol, "/*"• The definition 
of a data entry class may be restarted by the symbol 

In Step 3, the data set identification is requested. 
This entry normally requires identification of the infor- 
mation to be entered, its source and degree of processing, 
and appropriate reference coordinates. In the example, 
the test identification, "***TAM DEMONSTRATION***” , is 
entered. This step is terminated by an empty line. Step 4 
calls for the entry of attribute descriptions to further 
define quantitative values utilized in data entry. As an 
example, the entry ”A: Demonstration = 15", indicates 
that the attribute "A” has a value of 15 for the demonstra- 
tion. These definitions are terminated by an empty line. 

Step 5 begins the entry of quantitative data into 
the information system and the system calls for data entry 
definition. The first piece of information requested is 
the data entry type which is either point (P) , line (L) , 
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or area (A) information. In the example, point information 
is defined. The name of the data entry is requested next 
and a name "MMMM 11 is entered. Attribute names and values 
are requested from the operator with the attribute name 
first, and then the attribute value next. The attribute 
"A" is given a value of 15 in the example and attribute 
definition is terminated by the entry of an empty line. 

Next is requested the definition of data coordinates . At 
the conclusion of coordinate definition, an empty line is 
entered, and the operator is questioned as to whether he 
wishes to continue definition of additional data entries 
[Step 6). The response "N" terminates the define phase 
and completes the defined functions, whereas "Y" continues 
the process. 

Data Reconstruction 

An example of system operation is presented to 
demonstrate the handling of information data file point, 
line and area type of data. In Figure VII- 3 is shown a 
traditional insignia of Texas A§M University and the 
coordinates of the lines constructing the letters are 
tabulated. These coordinates are used in the system 
example and the letter T is given a value of 15, A has a 
value of 12 and M has a value of 4. In Figure VII-4 is 
shown the operator input sequence to input, transform, 
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Figure VII- 4. Operator Input Sequence for 
Computer Control 
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expand, and display these data on a computer generated 
printer map. 

The system directive XDATA allows operator specifica- 
tion of transformation of input data points. The output 
variable "<{>0 M is displayed, therefore the first attribute 
"A" is equated to "(j>0". In the example, the entry type, 
which corresponds to "AO" , is overridden and the END 
statement terminates the transformation definition. Input 
data entry types are overridden to produce point, line 
and area representations of the three letters. These 
printer maps are shown in Figures VII- 5, 6, and 7. 

Additionally, the system IRIMAS provides for coor- 
dinate translation, rotation and scaling through the 
XCOORD directive. In Figure VII- 8 is shown the operator 
definition sequence for the translation and 45-degree 
rotation of the university insignia. The result of these 
operations is shown in Figure VII-9. 

Brazos County Highway System 

The major highways in Brazos County, Texas represent 
a transportation resource which might be used in an anal- 
ysis of that county. In Figure VII- 10 is shown a general 
highway map of Brazos County. Highway systems represented 
on Hwy. 6, Hwy. 6 business, Hwy. 21, Hwy. 30, Hwy. 60, and 
Hwy . 1179 are entered into an IRIMAS data file, and 

^ ^ 
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Figure VII- 5. TAM Point Representation 
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Figure VII- 7. TAM 
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Figure VII- 8. Operator Definition for Rotation 
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Figure VII-9. TAM 45-Degree Rotation 
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represented by the attribute values 1-6, respectively. 

These highways are entered according to coordinates which 
represent the end point of straight line segments which 
appropriate the individual highways. A reconstruction of 
the highway system by IRIMAS from the piecewise linear 
segments is shown in Figure VII- 11. The application of 
data transformations provides for the selection and display 
of individual highway systems from the input files. In a 
like manner, coordinate transformations are applied to 
modify the scale or the orientation of the IRIMAS recon- 
struction to align the highway information with resource 
data from other sources . 

This demonstration of INFORM shows the construction 
of IRIMAS data files and the transformation and reconstruc- 
tion of data from these files. 

Another fruitful source of resource information is 
the imagery produced by the Earth Resources Technology 
Satellite. The next section addresses the creation of 
IRIMAS data files from interactive, IRIMAS- aided inter- 
pretation of ERTS imagery from Brazos County. 



Image Processing 

The software entitled IMAGE is that part of IRIMAS 
which implements the image analysis functions. This 
program segment contains the drivers for input and display 
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of ERTS data from digital magnetic tape. The image 
processing functions described in Chapter V are also 
implemented. Essentially, the same procedure required for 
the operator specification of function directives in the 
INFORM program segment is also required by IMAGE. 

The demonstration for the software segment IMAGE 
utilizes an ERTS scene from Brazos County, Texas. The 
portion of this scene to be analyzed is shown in Figure 
VII-12. This scene is of a region approximately 15 miles 
southwest of the city of College Station, Texas. The 
demonstration is devised to show typical procedures which 
might be utilized for the identification of the water 
resource represented by the Brazos River, and the incor- 
poration of this information into IRIMAS data files. 
Multiband imagery from ERTS is examined a single band at 
a time, classification parameters are determined through 
a training process, and the scene is classified to enhance 
the resource. The interpreted data is entered into an . 
IRIMAS data file for future reference. In Figure VII- 13 
is shown the CRT display system from which demonstration 
results are photographed. Radiance values measured by the 
ERTS satellite are stored on the digital magnetic tapes 
supplied by NASA, and these values are displayed in an 
image format on the CRT system. 
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Figure VII-12. BRTS Scene from Brazos County, Texas 
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Figure VII- 13. CRT Display System 
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IMAGE Operator Directives 

The operator control sequence for ERTS image data 
transformation and display is shown in Figure VII-14. 
Several features of IRIMAS are demonstrated here. The 
input driver for ERTS requests that the operator supply 

the computer tape line and cell starting coordinates and 
select the image band(s) to be analyzed. After the tape 
is positioned by the system at the start line, the program 
execution continues with a request for the definition of 
the data, transformation. The SLICE directive requires 
the system to prepare histogram and density slice infor- 
mation while the HISTOGRAM directive calls for a printer 
display of the data histogram which is also shown in 
Figure VII-14. 

Image Analysis 

Analysis of the ERTS scene typically begins with the 
preparation of data histograms for each of the ERTS image 
bands. This analysis step allows the operator to evaluate 
the range and distribution of the data values. Histograms 
consequently are .prepared for the Brazos Valley scene to 
be analyzed and are shown in Figure VII-15. The histogram 
printout gives the total number of points considered, and 
the number of picture elements at each radiance level. 

This information is also graphically displayed. 
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IRIMRS — RT4 INTERACTIVE RESOURCE 
INFORMATION AND ANALYSIS SYSTEM 
ENTER ANALYSIS FUNCTIONS 

0001 ERTS 

0002 TRANSFORM 

0003 SLICE 

0004 DISPLAY 

0005 60 TO 1 

0006 HISTOGRAM 
000? PAUSE 

ERTS-1 DATA INPUT DEFINE START LINE-' CELL 
LINES <0-276 0> > CELLS <0-81 0> 

START LINE* 825 
START CELL* 0 

SELECT BANDS < 1 « £ * 3 * 4 * A* ALL> 

BAND SELECT* A 

TAPE LINE NO. 0825 

ENTER .TRANSFORM 

00 01 I F < A4 - EQ . 0> THENA4* 128 

0 0 02 D 0»A 1 +A2 + 1 28-R3+ 1 28/ A4 

0003 A 1*00 16 

0004 END 

END OF^-DATA SET OH INPUT UN I T ♦♦CEJNT I NUF 

HI ST06RAM 
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16128 
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06220 
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Figure VII- 14. Operator Control Sequence 

for ERTS Processing 
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Histograms of ERTS Image Bands 
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Figure VII-15. 
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In Figure VI I- 16 are shown the CRT image displays of 
bands 1-3, which have been manually enhanced through inter- 
action with the CRT display to separate the radiance levels 
associated with the river. Band 4 is omitted because of 
the large amount of 0- level data indicated on the histogram 
which is attributed to poor data quality. 

Maximum likelihood classification provides for the 
resolution of elements of the scene which correspond to a 
defined category. In Figure VII- 17 is shown the operator 
sequence required for the training of the classifier 
software. Upon execution, the system requests the number 
of classes to be considered, and the coordinate boundaries 
of the rectangular areas from which the data are to be 
considered. Names of the data sets are additionally 
requested. The STATISTICS directive causes the printouts 
of the results of the training phase. These results for 
the example are shown in Figure VII- 18. The training 
phase prepares the statistical parameters required by the 
classification algorithm. 

The operator directive sequence for the specifica- 
tion of maximum likelihood classification of ERTS imagery 
is demonstrated in Figure VII- 19. The classifier portion 
of IRIMAS utilizes the class statistics obtained during 
the training phase for classification. The result of this 
classification with only the "water" class displayed is 
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Figure VI I- 16. Display of ERTS Image Bands 




202 


» 



Figure VI I- 16. Display of ERTS Image Bands (cont.) 
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ENTER ANALYSIS FUNCTIONS 

0001 ERTS 

0002 TRAIN 

0003 60 TO 1 

0004 PAUSE 

ER7S-1 DATA INPUT DEFINE START LINE^CELL 
LIMES < 0-2760} ? CELLS <0-31 0> 

START LINE* 325 
START CELL* 0 

SELECT BANDS <1 * £ * 3 ? 4 > A=ALL > 

BAND SELECT* A 
TAPE LINE NO. 0325 
CLASSIFICATION TRAINING SETS 
ENTER NUMBER OF CLASSES <i-4> 4 


NAME OF CLASS CIS W 1 

LEFT-HAND BOUNDARY: 32 

RIGHT-HAND BOUNDARY: 93 

UPPER BOUNDARY: 94 

LOWER BOUNDARY: 96 

NAME OF CLASS C£: W2 

LEFT-HAND BOUNDARY: 122 

RIGHT-HAND BOUNDARY: 127 

UPPER BOUNDARY: 126 

LOWER BOUNDARY: 134 



NAME DF CLASS C3: 
LEFT-HAND BOUNDARY: 
RIGHT-HAND BOUNDARY: 
UPPER BOUNDARY: 

LOWER BOUNDARY: 


W3 

133 

143 

133 

141 


NAME OF CLASS C4: W4 

LEFT-HAND BOUNDARY: 140 

RIGHT-HAND BOUNDARY: 157 

UPPER BOUNDARY: 174 

LOWER BOUNDARY: 183 


END OF DATA SET ON INPUT UNIT^CONTINUE 


Figure VII- 17. Operator Sequence for Classifier Training 
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CLASS 1 
TOTAL POINTS 

0 0 03 3 

MEAN VALUES 

0032 0029 0043 OOZ'D 
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figure VI I- 18. Classifier Training Results 
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ENTER ANALYSIS FUNCTIONS 
0001 ERTS 
000c CLASSIFY 

0003 TRANSFORM 

0004 DISPLAY 

0005 GO TO 1 

0006 PAUSE 

ERTS-1 DATA INPUT DEFINE START LINE-CELL 
LINES < 0-376 0 > ? CELLS < 0-3 1 0> 

START LINE- 825 
START CELL- 0 

SELECT BAND'S \ 1 j 2 » 3 ? 4 ? ft—ALL - 1 
BAND SELECT- 2 
TAPE LINE NO. 0825 
ENTER TRANSFORM 
0001 IF<fll.NE.. 

8> THEMft 1 = 0 

0 0 02 I F < A 1 . EQ . 2> THENA 1 = 15 
0003 END 

END DF DATA SET ON INPUT UN IT ♦♦CONTI HUE 
DEFINED FUNCTIONS APE COMPLETE ©♦♦CON 7. 
ENTER ANALYSIS FUNCTIONS 

0001 



Figure VII- 19. Operator Specification of 

Classify Sequence 
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shown in Figure VII- 20. At this point in the analysis, 
the water resource identified as the Brazos River through 
interpretation of the image presentation is interactively 
saved in the IRIMAS data files. In Figure VII- 21 is shown 
the listing of this data file. 

Change Monitoring 

A third demonstration is used to illustrate the 
versatility of the IRIMAS data system in the monitoring 
of resource changes over a period of time. Monitoring of 
changes using the traditional grid format storage would be 
exceedingly difficult especially with the overlay and 
registration of mult i- temporal images. Use of the IRIMAS 
data formats and analysis capabilities make this monitoring 
a simple interactive task. 

The feature chosen to show a procedure for change 
analysis is the west-bypass of U.S. Highway 6 around the 
Bryan- College Station, Texas communities. Earth Resources 
Technology Satellite data from the area has been accumu- 
lated for the period August 1972 through April 1973, a 
period during which the bypass has shown considerable 
development. The latest available general highway map of 
the area (see Figure VII-10, page 190) does not show the 
west bypass, however, construction appears on the recent 
ERTS imagery. In Figure VII- 22, are shown several image 
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Figure VII-20. Classification Results for "Water" Class 
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Figure VII-21. IRIMAS File Listing of Brazos River Data 
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Figure VII- 22. ERTS Tmage Sequence During Highway 6 

West Bypass Construction 
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Figure VII- 22. ERTS Image Sequence During Highway 6 

West Bypass Construction (cont.) 
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reproductions of the area under consideration for the 
period where development of the bypass is apparent. Bypass 
changes are interpreted from the imagery and are added to 
the IRIMAS data file created for the general Brazos County 
highway system. An IRIMAS expansion of the data file is 
displayed in Figure VII-23. 

The sequence of highway construction is indicated by 
different letters on the printout. The symbol 
represents the major highway systems which occur on the 
general highway map. Bypass changes interpreted from 
ERTS imagery are encoded by date as follows: 1) B- 5/27/73, 
2) C- 7/20/73, 3) G-7/2G/73, 4) D- 3/29/74. 

Maintenance of the interpreted information requires 
minimal effort in IRIMAS data files where changes in 
resource information can be readily saved. Reconstruction 
of these data is also provided for by the IRIMAS system. 
Traditional data formats and image interpretation systems 
preclude the simple maintenance of the interpreted data 
in a single file. Full images must be maintained and 
overlayed for each date considered. The overlay method 
also requires full interpretation of the image for each 
analysis and does not take advantage of previous inter- 
pretations, whereas IRIMAS provides for accurate recall 
of past interpretations. 
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CHAPTER VII 

CONCLUSIONS AND RECOMMENDATIONS 

The recent availability of timely synoptic earth 
imagery from the Earth Resources Technology Satellite 
(ERTS) provides a wealth of information for the monitoring 
and management of vital natural resources. Satellite data, 
as well as information derived from aircraft sensors, are 
providing an increasingly more important supplement to 
conventional data sources used for the monitoring, protec- 
tion, and allocation of resources. Remote sensor data 
provides a unique point of view which enhances the infor- 
mation obtained about an area resource from traditional 
ground survey sources. It is becoming apparent, due to 
the scope and amount of data available, that computer 
processing and analysis of resource information will 
become a necessity for timely application. Computer soft- 
ware systems have not existed to provide the planning 
function with an adequate information structure to main- 
tain and utilize spatial and temporal data interpretations. 
Available systems, for the most part, also lack the sim- 
plistic interface required by the resource planner who is 
not versed in the complexities of computer processing. 

The resource management decision is concerned with 
the evaluation of extrapolated land usage and its 
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associated environmental impact, both in an assumed steady 
state and in the presence of actual or anticipated pertur- 
bations. The resource manager asks composite questions 
about the relationships of area extensive resource fea- 
tures. These questions amount to a conceptual overlay of 
any of the ’'maps” which comprise the data base. Decisions 
are based on the disposition of related resources such as 
soils, water, and vegetation. 

The objective of this research has been to develop 
an information management system suited to the general 
needs of resource managers engaged in the utilization of 
imagery and other spatiotemporal data in their management 
processes. The need for such a system has been inferred 
by numerous authors [1], [2], [6], [26] who have empha- 
sized the need for an interactive analysis system capable 
of handling the information requirements of the spatial 
and temporal data available from a variety of sources. 

Several processing systems have been developed to 
address one aspect or another of the resource management 
problem, however each of these fails with respect to one 
or more of the essential ingredients for an effective 
management tool. These ingredients include a simplistic, 
interactive interface; the flexibility to handle a variety 
of information types; and a system whose costs are in-line 
with operating budgets at a regional agency level. The 



215 


Interactive Resource Information Management and Analysis 
System (IRIMAS) has been developed to specifically 
address the need for each of these ingredients in the 
handling and analysis of the spatial and temporal informa- 
tion required for resource related decision-making. The 
system provides for the interactive interpretation of 
data from ERTS imagery, cartographic products, and point 
surveys. Operator interface is simplistically accomplished 
through single word commands and computer initiated queries 
and requires only a minor amount of operator training for 
its effective use. A flexible data base has been set up 
to address point, line and area format data which can be 
retrieved for the display of overlays according to user 
priority-weights. Additionally, the system can be dupli- 
cated for equipment costs below $50,000, or may be set up 
in an interactive time-shared mode to take full advantage 
of existing machines. 

IRIMAS is innovative from several significant aspects. 
The application of formal language theory to the manage- 
ment of information data files has been alluded to only 
recently in the literature, yet it provides the basis for 
the flexible information format developed for this system. 
The transformation and maintenance of interpreted sensor 
information in data files, compressed by the utilization 
of category boundaries, offers a worthwhile alternative to 
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the accepted processing standards prevalent in image 
analysis today. The greatest significance of IRIMAS, 
however, lies in the evolution of a new approach to the 
utilization of remote sensor data and computer processing 
techniques in resource analyses and in the introduction of 
a new tool which exhibits a great amount of flexibility 
and potential for interactive interpretation of resource 
information. 

The Interactive Resource Information Management and 
Analysis System, IRIMAS, is a hardware/software system 
developed with a Texas Instruments 980A minicomputer sys- 
tem consisting of a CPU with 12K words of memory, a 9-track 
digital magnetic tape unit, a send-receive terminal with 
dual magnetic tape digital cassettes, and a color CRT 
image display. Several additional pieces of equipment 
would suitably augment this hardware arrangement. Exper- 
ience shows that the 12K words of memory is inadequate to 
handle the fully implemented IRIMAS software package. It 
is recommended that future configurations of IRIMAS can be 
provided with an additional 8K words of memory to add 
system flexibility, and to remove some of the software 
restrictions placed on the system. 

The information storage format requires the defini- 
tion of resource features in terms of coordinate bounda- 
ries. It is time-consuming and cumbersome to manually 
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extract the coordinate information and to enter it into 
the data files. This defining operation could be more 
efficiently accomplished with the aid of an automated 
graphics input device. Several varieties of these instru- 
ments are available commercially and would greatly enhance 
the efficiency of IRIMAS utilization. Future implementa- 
tion of IRIMAS could effectively utilize this class of 
input device, and its inclusion in the system is rec- 
ommended. 

Several inadequacies are noted in the IRIMAS software 
and modifications are recommended. The cluster algorithm 
implemented in IRIMAS has produced less than desirable 
results in the analysis of ERTS imagery.. Several reasons 
are suggested for this lack of performance. Primarily, 
the difficulty seems to lie in the number of classes 
permitted, which has been limited to sixteen, because of 
memory constraints. Imagery from ERTS characteristically 
has more individual classes associated with a scan line 
than was experienced with data sets in previous applica- 
tions of the algorithm. These data differences suggest 
that cluster thresholds should also be examined for the 
ERTS data. 

The maximum likelihood classification algorithm 
applied in IRIMAS, has a fixed threshold for the identifi- 
cation of the "no classification” category which has been 
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set high enough to allow classification of all points. It 
is recommended that this threshold be set by the system 
operator with the specification of percentage of points to 
be included or rejected by the algorithm. Computation or 
table look-up could be applied for the determination of 
the appropriate threshold as a function of the number of 
dimensions and individual class statistical parameters. 

This procedure would facilitate the single category class- 
ification of a scene with the rejection of other categories 
which would be useful in a resource interpretation system. 

Coordinate and data transformation software is imple- 
mented to support fixed-point binary arithmetic. This 
method was implemented to take advantage of the hardware 
integer multiply/divide feature of the TI980A system. 
Operational use of the system has indicated that the 
transformations normally required would be more flexible 
and precise if arithmetic operations were performed by 
floating-point operations rather than the fixed-point 
binary operations. It is recommended therefore that the 
transformation interpretor be modified to perform floating- 
point arithmetic rather than fixed-point arithmetic. 

Operational use of the IRIMAS classification algo- 
rithms is extremely slow computationally. Classification 
of the 128 x 128 pixel scene using a single ERTS band 
requires in excess of 3 minutes, and use of the full four 
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bands expands the time to 12-15 minutes. These times are 
undesirable from an interactive point of view and could be 
significantly reduced if the classification software were 
converted to assembly language rather than FORTRAN. 

Other system improvements could be accomplished by the 
incorporation of software drivers to allow inclusion of 
additional image format data sources, and to include image 
spatial- frequency processing. With the addition of a disc 
to the TI980A system, drivers could be developed to use the 
sector addressable capability of that system to enlarge the 
processing window by writing to the disc rather than into 
computer memory. Overlay link-editing and real-time load- 
ing procedures could be applied in the disc system to im- 
prove the efficiency of memory usage. Additionally, it 
seems practical to develop the software to allow utiliza- 
tion of the CRT graphics-overlay feature to superimpose 
resource information graphically into an image display. 
These recommended systems improvements should add consider- 
able versatility to the system and should enable more 
effective use of the system as an interactive resource 
management aid. 
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